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Two-dimensional (2D) transition metal dichalcogenide (TMD) nanosheets exhibit remarkable elec¬ 
tronic and optical properties. The 2D features, sizable bandgaps, and recent advances in the syn¬ 
thesis, characterization, and device fabrication of the representative M 0 S 2 , WS 2 , WSe 2 , and MoSe 2 
TMDs make TMDs very attractive in nanoelectronics and optoelectronics. Similar to graphite 
and graphene, the atoms within each layer in 2D TMDs are joined together by covalent bonds, 
while van der Waals interactions keep the layers together. This makes the physical and chemical 
properties of 2D TMDs layer dependent. In this review, we discuss the basic lattice vibrations of 
monolayer, multilayer, and bulk TMDs, including high-frequency optical phonons, interlayer shear 
and layer breathing phonons, the Raman selection rule, layer-number evolution of phonons, mul¬ 
tiple phonon replica, and phonons at the edge of the Brillouin zone. The extensive capabilities of 
Raman spectroscopy in investigating the properties of TMDs are discussed, such as interlayer cou¬ 
pling, spin-orbit splitting, and external perturbations. The interlayer vibrational modes are used 
in rapid and substrate-free characterization of the layer number of multilayer TMDs and in probing 
interface coupling in TMD heterostructures. The success of Raman spectroscopy in investigating 
TMD nanosheets paves the way for experiments on other 2D crystals and related van der Waals 
heterostructures. 


I. INTRODUCTION 

Driven by the unique properties of graphene, the rapid 
progress of graphene research has stimulated interest in 
other two-dimensional ( 2 D) crystals.®® Several represen¬ 
tative 2 D crystals, such as BN, M0S2, NbSe2, B^Tes, 
Bi 2 Se3, and B^S^CaC^O^,® have been obtained by 
similar methods to those used to obtain graphene. The 
library of 2 D crystals is very rich, including many natu¬ 
rally occurring layered materials (LMs) (e.g., minerals), 
transition metal oxides, and transition metal dichalco¬ 
genides (TMDs) .® The capability to prepare/grow differ¬ 
ent families of 2 D crystals will greatly enhance the library 
of materials, and thus allow exploration of the fascinating 
unusual physics in two dimensions.® Recently, 2 D TMDs 
such as M0S2, MoSe2, WS2, and WSe2 have attracted 
considerable attention because of their unique properties 
and their potential applications in (opto)electronics PflU 
These 2 D crystals can be grown or mechanically ex¬ 
foliated to monolayer thickness, similar to the exfolia¬ 
tion of graphene. However, in contrast to graphite and 
graphene, TMDs consist of more than one element (as 
shown in the case of M0S2 in Fig. [lja), where one 
Mo plane is sandwiched between two S planes), which 
makes their lattice dynamics more complex than mul¬ 
tilayer (ML) graphene ,® 7 including the symmetry, force 
constants, and frequency variation with thickness EBHUI 
Furthermore, compared with the bulk form, monolayer 
( 1 L) and few-layer (FL) TMDs show very distinctive 
physical properties, especially the lattice dynamics, elec¬ 
tronic structure, and optical properties. 

Many TMDs possess a graphene-like layered structure 
with the single carbon atomic layer of graphene replaced 
by the repeat unit X-M-X, where M is a transition metal 
of group VI (Mo or W) and X is a chalcogen (S, Se, or 


Te). More than 40 different types of TMDs can occur 
depending on the combination of chalcogen and metal. 
For example, in M0S2, the chalcogen atoms (S) in the 
two hexagonal planes are separated by a plane of metal 
atoms (Mo), as shown in the case of IL-M0S2 in Fig. 
1 ^a). As indicated by the black triangular prism in Fig. 
l'a), the Mo atom has trigonal prismatic coordination, 
where the upper triangle is obtained by mirroring the 
bottom triangle. This repeat unit stacking along the c 
axis will form two different polytypes, i.e., 2 H (hexago¬ 
nal symmetry, two X-M-X layers per repeat unit, trigo¬ 
nal prismatic coordination, H stands for hexagonal) and 
3 R (rhombohedral symmetry, three X-M-X layers per re¬ 
peat unit, trigonal prismatic coordination, R stands for 
rhombohedral).® The site of the metal atom can also have 
octahedral coordination, where the upper triangle is in¬ 
version of the_bottom triangle, which gives the third poly¬ 
type 1 T F I 22 I 23 I (tetragonal symmetry, one X-M-X unit in 
the unit cell, octahedral coordination, T stands for trig¬ 
onal). Schematic diagrams of the three structure poly¬ 
types are shown in Fig. [ljb).® The difference between 
the trigonal prismatic symmetry and octahedral coordi¬ 
nation of metal atoms can result in completely differ¬ 
ent electronic structures in the monolayer. For example, 
1L 2H-M0S2 is semiconducting,®^ while 1L IT-M0S2 is 
metallic P 3 * 24 ! Recently, 1L IT-M0S2 obtained by interca¬ 
lation with alkali metals transformed to 2H-M0S2 by an¬ 
nealing at 300 °C, which was verified by the re-emergence 
of bandgap photoluminescence (PL)J^ The symmetry 
difference between the monolayer of 2H-MX2/3R-MX2 
(point group D^h) and IT-MX2 (point group D^d) can 
be determined by Raman analysis.^ Here, we mainly fo¬ 
cus on 2 H-MX 2 , and bulk 2 H-MX 2 is simply denoted as 
bulk MX 2 . 

Great advances have been made with M0S2, and re- 
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FIG. 1: (a) Schematic diagram of the three-dimensional structure of IL-M 0 S 2 . The triangular prism (black lines) shows the 
trigonal prismatic coordination of Mo atoms, (b) Schematic diagrams of the three typical structural polytypes of MX 2 : 2H, 3R, 
and IT. a and c represent the lattice constants. Reproduced with permission from ref. 0 ] Copyright 2012, Nature Publishing 
Group, (c) Optical image of 1L- and 2L-MoS 2 obtained by micromechanical cleavage, (d) Stable suspensions of MoS 2 , WS 2 , 
and BN from liquid-phase exfoliation in solvents (top), and thin films derived from vacuum filtration of these suspensions 
(bottom). Reproduced with permission from ref. 1211 Copyright 2013, American Association for the Advancement of Science, 
(e) Scanning electron microscopy images of triangular MoS 2 thin films by vapor phase growth. Reproduced with permission 
from ref. 12T1 Copyright 2013, Nature Publishing Group. 


search is gradually being extended to other TMDs. Most 
TMDs are indirect-gap semiconductors in the bulk, but 
transform to direct-ga p sem iconductors when reduced to 
monolayer thickness! 8 * 9 * 26 * This can be probed by a va¬ 
riety of spectroscopic tools.® 9 ^HS] For example, as an 
indirect-gap material, bandgap PL in bulk 2H-M0S2 is 
very weak because it is a phonon-assisted process and 
known to have negligible quantum yield. Appreciable 
PL is observed in FL-M0S2 and surprisingly bright PL 
is detected in IL-M0S2, which is indicative of it being a 
direct-gap semiconductor.®^ The indirect-to-direct tran¬ 
sition is attributed to the absence of weak interlayer cou¬ 
pling in the monolayer. 8 9 ^* which was directly verified 
in FL-MoSe2 via interlayer thermal expansion. 33 The di¬ 
rect bandgap of monolayer TMDs promises high perfor¬ 
mance in optoelectronic applications, such as the extra- 
high on/off ratio reported for M0S2 transistors at room 
temperature .^ 0 Two peaks (A and B) are observed in the 
absorption spectra of 1L- and 2L-M0S2, corresponding to 
transitions between split valence bands (VBs) and con¬ 


duction bands (CBs)P The energy splitting between A 
and B in bulk M0S2 arises from the combined effect of 
interlayer coupling and spin-orbit coupling (SOC), but 
solely comes from SOC in IL-M0S2 because of the ab¬ 
sence of interlayer couplingP Furthermore, because of 
the strong SOC and inversion symmetry breaking in 1L- 
M0S2, coupled spin and valley physics introduces a new 
controllable degree of freedom beyond charge and spin, 
leading to valleytronicsP^X 16 An excited electron will at¬ 
tract a hole by Coulomb interactions to form a bound 
state (called an exciton) with the bound energy levels 
lying in the bandgap region. This excitonic effect is 
significantly enhanced in low-dimensional materials ow¬ 
ing to stron g sp atial confinement and a reduced screen¬ 
ing effect P™^* Thus, as a 2D material, 1 L-MX 2 is ex¬ 
pected to have a strong excitonic effect with a large ex¬ 
citon binding energy of 0 . 5 - 1.0 eV®®* which will signif¬ 
icantly influence its optical properties. In addition, one¬ 
dimensional nonlinear optical edge states are observed at 
the edge of IL-M0S2, which allows direct optical imaging 
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of the atomic edges and boundaries of the 2 D material.® 

Various methods have been developed to obtain uni¬ 
form, high-purity, thickness-controlled, and large-area 
TMDs monolayers, which can be mainly divided into 
two categories: exfoliation from parent bulk crystals and 
chemical vapor deposition (CVD) on substrates. The 
coupling between the X-M-X repeat units of bulk TMDs 
by weak van der Waals forces allows the crystal to be 
easily cleaved along the layer surface. Methods inherited 
from th e grap hene community, such as micromechanical 
cleavage j 18 l 45 > liquid-phase exfoliation, 3 ® and exfoliation 
by intercalation with ionic species,® have been widely 
used in preparing monolayer TMDs. Samples obtained 
by micromechanical cleavage exhibit high crystal quality, 
but need to be supported on substrates (e.g., Si 0 2 /Si, 
as shown in Fig. me)) for optical identification. 4 -^® 
Liquid-phase preparation of MX 2 (Fig. [ljd)) is promis¬ 
ing, and allows the pr epar ation of composites and hy¬ 
brids of different MX 2 P® However, the thickness and 
size of thin layers cannot be well controlled compared 
with exfoliation processes. Inspired by CVD of large-area 
graphene on copper,® several CVD methods based on in¬ 
sulating substrates have recently been reported!®® 52 ^ 
CVD growth of TMD monolayers offers the possibility of 
wafer-scale fabrication of electronic and optoelectronic 
devices!® 

Ideally, methods for fundamental research and ma¬ 
terial characterization should be fast and nondestruc¬ 
tive, offer high spectral and spatial resolution, provide 
structural and electronic information, and be applica¬ 
ble at both laboratory and mass-production scales. Ra¬ 
man spectroscopy fulfills all of the above requirements. 
Raman spectroscopy can identify unwanted byproducts, 
functional groups, structural damage, and chemical mod¬ 
ifications introduced during preparation, processing, or 
placement of various LMs. Here, we review recent de¬ 
velopments in Raman spectroscopy for characterization 
of 2D TMDs. We first introduce lattice vibrations and 
phonon dispersion of 1 L, 2 L, and bulk MX 2 (Section 2 ). 
Then, we discuss the fundamentals of Raman scattering 
and its selection rule (Section 3 ). We outline nonreso¬ 
nant (Section 4 ) and resonant Raman (Section 5 ) spec¬ 
troscopy of TMDs. In Section 6, Raman spectroscopy 
of 2D TMDs modified by external conditions, such as 
strain, pressure, temperature, and electric field, is dis¬ 
cussed. Finally, we discuss the potential of shear and 
layer breathing modes for characterization of the thick¬ 
ness of unknown ML TMDs and probing the interface 
coupling of 2D hybrids and heterostructures (Section 7 ). 
We hope that this review will enable a better under¬ 
standing of the recent progress in phonon and Raman 
scattering of 1 L, ML, and bulk TMDs, and be helpful to 
design and perform Raman experiments to make use of 
the great potential of Raman spectroscopy for LMs, and 
their hybrids and heterostructures. 


II. LATTICE VIBRATIONS AND PHONON 
DISPERSION OF 1L-, 2L- AND BULK 2H-MX 2 

A. Lattice vibrations at the T point and its 
symmetry classification 

Lattice vibrations can be classified based on the ir¬ 
reducible representation of the symmetry group of the 
crystals.® The unit cell of bulk MX 2 consists of two X- 
M-X units with a total of six atoms, suggesting that there 
are 18 phonon modes (3 acoustic and 15 optical modes). 
Bulk MX 2 has D 6 h point group symmetry. Thus, the 
lattice vibrations of bulk MX 2 at T can be ex presse d by 
the irreducible representations of Dq^ as follows— 57 : Y— 
A\ g -\-2A 2u -\-2B 2g -\-Bi u -\-E\ g -\-2Ei u -\-2E 2g -\-E 2u , where 
one A 2u and one E\ u are acoustic modes, A\ g , Ei g , and 
E 2g are Raman (R) active, another A 2u and E\ u are in¬ 
frared (IR) active, and B 2g , Bi ui and E 2u are optically 
inactive (silent). Here, the modes denoted by the letter 
n E” are doubly degenerate in the xy plane. 

Because of the lack of translational symmetry (r) along 
the z axis, which is perpendicular to the basal xy plane, 
there is a reduction in symmetry in FL TMDs (NL- 
MX 2 ),® where N is an integer and refers to the layer 
number. Odd NL-MX 2 (ONL-MX 2 ) has D 3 ^ point group 
symmetry because of the presence of the horizontal reflec¬ 
tion plane (07J that passes through the transition metal 
atom (M). In particular, the unit cell of 1 L-MX 2 is com¬ 
posed of three atoms with nine normal vibrational modes 
at the T point, which is expressed based on D3^ 

T= 2A 2 +A 1 +2E +E , where one A 2 and one E are 
acoustic modes, another A 2 is IR active, A 1 and E are 
R active, and another E is both R and IR active, as 
shown in Fig. [ 2 |a). The point group of even NL-MX 2 
(ENL-MX 2 ) is D3^ because of the existence of inversion 
symmetry i. ah exists in ONL-MX 2 and bulk MX 2 but 
it is not present in ENL-MX 2 .® In particular, for 2 L- 
MX 2 , the unit cell is composed of six atoms belonging 
to D 3 d. There are 18 normal vibrational modes at the 
r point!- 20 ® T= 3Aig+3A 2u +3E g +3E u , where one A 2u 
and one E u are acoustic modes, the other A 2u and E u 
are IR active, and A\ g and E g are R active, as shown in 
Fig. [ 2 ja). It should be noted that the R and IR modes 
are mutually exclusive in MX 2 and ENL-MX 2 because 
of the presence of inversion symmetry i. The symme¬ 
try and normal mode displacement of each mode of bulk, 
2 L, and 1 L MX 2 are shown in Fig. Ufa). In general, the 
9 N normal vibrational modes in NL-MX 2 at the Y point 
can be expressed based on its corresponding point group 
a j2U | 58 | 59 | YoNL = (A[ + E" ) + + e') for 

ONL-MX2 (D3/j) and Yenl = ^f(Aig+E g + A2 U + E u ) 
for ENL-MX 2 (D 3 d). These 9 N modes can be further 
classified into three categories according to their fre¬ 
quencies (see Section 4 ): acoustic, ultralow frequency 
(ULF), and high frequency modes. These classifications 
are shown in Table. [l| where the number of Raman active 
modes for the three categories are also included. Like- 
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FIG. 2: (a) Symmetry and normal displacements of each optical vibration mode for bulk MX 2 , 2 L-MX 2 , and IL-MX 2 . The 

Raman-active (R), infrared-active (IR), and both R and IR inactive (silent) modes are identified. The dotted lines show that 
each mode in IL-MX 2 is split into two modes in both 2 L-MX 2 and bulk MX 2 , of which one vibrates in-phase and another out- 
of-phase. Calculated phonon dispersion curves of (b) IL-MX 2 and (c) bulk MX 2 for M 0 S 2 . The points in (c) are experimental 
data extracted from the neutron scattering measurements in ref. 1531 Reproduced with permission from ref 1541 Copyright 
2011, American Physical Society. 


wise, the symmetry classifications for lattice vibrations 
away from the T point, e.g., at high-symmetry edges of 
the Brillouin zone, can also be obtained according to the 
corresponding point group, which has lower symmetry 
than that at the T point 

Each of the nine normal vibrational modes in IL-MX 2 
will split into the corresponding two modes in both 2L 
and bulk MX 2 , as shown by the dotted lines in Fig. 
[ija). For example, E' (R+IR) in 1L-MX 2 splits into 
Eg (R) and E u (IR) in 2L-MX 2 , or E\ g (R) and E\ u 
(IR) in bulk MX2, where the displacement between the 
top and bottom layers are in phase for E u /Ei u: but out 
of phase for E g / E 2g . The positive integer in the upper 
right corner of the mode notation is used to distinguish 
modes with the same symmetry (e.g., E\ g and E 2g in 
Fig. § a)), and is generally applicable to other modes. 


In the case of E g /E 2g in 2L and bulk MX 2 , additional 
’’spring” coupling between X atoms of neighboring lay¬ 
ers increases the frequency with reapect to E u /E\ u . This 
is known as ’’Davydov splitting”which also occurs to 
other modes in 2L and bulk MX 2 . Thus, the position of 
E\ g (Pos(£ 1 2 ^)) in bulk MX 2 is expected to be larger than 
the position of E\ u (Pos(Fi n )). However, this is not the 
case in bulk M 0 S 2 , where Pos^^) is ~1 cm -1 lower than 
Pos(£’i u ) because of the existence of causes other than 
the weak interlayer interaction. 61 Interesting splitting oc¬ 
curs in acoustic A 2 (ZA) and E (LA/TA) in 1L-MX 2 . 
A 2 (ZA) in IL-MX 2 will split into A\ g and A 2u (ZA) 
in 2L-MX 2 , or B 2g and A 2u (ZA) in bulk MX 2 , while 
E (LA/TA) in IL-MX 2 splits into E g and E u in 2L- 
MX 2 , or E\ g and E lu (LA/TA) in bulk MX 2 . Al g /B$ g 
and E g /E 2g in 2L and bulk MX 2 are referred as layer 
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TABLE I: Symmetry, Raman activity, and number of acoustic, ultralow frequency, and high frequency modes in ENL, ONL, and 
bulk MX 2 . R, IR and S in brackets after the symmetry refer to Raman active, infrared active, and optically silent, respectively. 
All of the E modes are doubly degenerate. Num(R) denotes the number of Raman active modes. 


NL-MX 2 

Acoustic 

Ultra-low frequency 

High frequency modes 

Num(R) 



C modes 

LB modes 



Bulk 

A2uAE\ u 

E 2 2g (R) 

B%{S) 

Eig(R) + E\ U (IR) + E 2g {R) + E 2 u {S) 
+A lg (R) + A 2u (IR) + B 2g (S) + Bi u (S) 

7 

ENL 

A 2 uT E u 

f E g (R)+EE> e u {IR) 

f A lg {R)+^l A 2u (IR) 

N(E g (R) + E U (IR) + A lg (R) + A 2u {IR)) 

9N/2 

ONL 

A 2 -\-E 

Zf±(E(IR+R) + E"(R)) 

V(^'i (R) + A 2 (IR)) 

N(E / (IR + R) + e"{R) + A 1 (R) + A 2 (IR)) 5(3N-l)/2 


breathing (LB) and shear modes, which corresponds to 
the relative motions of two X-M-X layers perpendicular 
and parallel to the layer plane, respectively. A similar 
shear mode has been observed in multilayer graphene. 17 
It is usually called C because it provides a direct mea¬ 
surement of the interlayer coupling. 17 


B. Phonon dispersion 

Away from the T point, normal vibrational modes are 
dispersive with respect to wavevector q. The phonon 
dispersions of 1L and bulk M 0 S 2 calculated by density 
functional theory within the local density approximation 
are shown in Fig. ^b) andj^c)!^ 

The phonon dispersion of IL-MX 2 has three acoustic 
and six optical branches inherited from the nine vibra¬ 
tional modes at the T point. The three acoustic branches 
are the in-plane longitudinal acoustic (LA), the trans¬ 
verse acoustic (TA), and the out-of-plane acoustic (ZA) 
modes. The LA and TA branches have linear disper¬ 
sion and a higher frequency than the ZA mode around 
T.E3E2 63 The six optical branches are two in-plane longi¬ 
tudinal optical (LOi and LO 2 ), two in-plane transverse 
optical (TOi and TO 2 ), and two out-of-plane optical 
(ZOi and ZO 2 ) branches. These six optical modes at 
the T point correspond to the irreducible representations 
(Fig. H(b)) E" (LOi and TOi), E (L0 2 and T0 2 ), 
A 2 (ZOi), and A 1 (ZO 2 ). Because MX 2 compounds are 
slightly polar materials, certain IR-active phonon modes 
(e.g., E ) show LO-TO splitting because of the coupling 
of the lattice to the macroscopic electric field created by 
the relative displ acement of M and X atoms in the long- 
wavelength limit IM1M1MJ i n addition, there is a bandgap 
between the acoustic and optical branches, which is ^100 
cm -1 for M 0 S 2 (Fig. [2^b) and[2jc)) and WS 2 , ^30 cm -1 
for WSe 2 , and ~15 cm -1 for MoSe 2 (see Fig. SI in the 
Electronic Supplementary Information (ESI)) ! 54 l 62 l 63 l 

For 2L and bulk MX 2 , there are 18 phonon branches, 
which are split from the nine phonon branches in 1L- 
MX 2 . Owing to the weak van der Waals interlayer in¬ 


teractions in 2L and bulk MX 2 , the frequency splitting 
of the corresponding two branches in 2L and bulk MX 2 
from each optical branch in IL-MX 2 is very small, result¬ 
ing in similar optical phonon dispersion curves for 1L, 2L, 
and bulk MX 2 . For example, the phonon dispersions of 
1L (Fig. |2jb)) and bulk M 0 S 2 (Fig. [2|c)) are very simi¬ 
lar, except for the three new branches below 100 cm -1 in 
bulk M 0 S 2 because of interlayer vibrations. The general 
features of 1L, 2L, and bulk MoS 2 P WS 2 ! ^§3 WSe 2 P 
and MoSe 2 -® are similar. However, their phonon frequen- 
cies are very different (see Fig. [2jb) and Fig. SI in ESI). 
Compared with IL-M 0 S 2 , the phonon bands in IL-WS 2 
are shifted to lower frequencies, which is mainly because 
of the larger mass of W atomsSi milar behavior oc¬ 
curs between lL-WSe 2 and IL-WS 2 ! 63 * 66 ! Furthermore, 
the interlayer phonons in bulk MX 2 are different, which 
is because of the intrinsic differences in both the complete 
mass of MX 2 and the interlayer coupling strengths. 


III. RAMAN SCATTERING AND ITS 
SELECTION RULE FOR PHONONS IN 1L-, FL- 
AND BULK MX 2 

Raman scattering is one of the main methods to obtain 
information about the lattice vibrations of a crystal. In 
the Raman scattering process for a crystal, an incident 
photon from a laser can absorb or emit a phonon of the 
crystal, and is then scattered into the surroundings in all 
directions. As a result, the photon loses (Stokes shift) or 
gains (anti-Stokes shift) the energy of the phonon, which 
is then detected and analyzed. Such a Raman process 
is called a first-order Raman process. The process is 
called a second-order Raman process if two phonons are 
involved. Under the restriction of momentum conserva¬ 
tion, only the phonon with wavevector q « 0 (T point) 
can contribute to first-order Raman scattering. The two 
phonons in the second-order Raman process should have 
opposite wavevector q. Second-order Raman scattering 
is usually much weaker than first-order scattering, but 
will be resonantly enhanced by the interband transition 
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FIG. 3: Calculated Raman spectra of (a) bulk M 0 S 2 and 
(b) IL-M 0 S 2 with polarization vectors of e*(l, 0, 0) and e s 
(cos #, sin #, 0) for the incoming and scattered photons, respec¬ 
tively. The inset figures show I(El g )/I(A lg ) and 1{E )/I{A 1 ) 
as a function of #. Reproduced with permission from ref. [671 
Copyright 2014, Royal Society of Chemistry, (c) Raman spec¬ 
tra of 5L, 6L, and bulk M 0 S 2 measured under XX (red) and 
YX (black) polarization configurations. Here, XX and YX 
refer to 0 = 0 and 6 — 7r/2, respectively. The irreducible 
representation of each mode is indicated. Reproduced with 
permission from ref. 1191 Copyright 2013, American Physical 
Society. 


of a semiconductor. This enhancement makes it possi¬ 
ble to detect the phonons at the edge of the Brillouin 
zone (e.g., M and K points).^ Furthermore, under the 
excitonic resonance condition, the symmetry of exciton 
levels can mediate the scattering of the Raman-inactive 
phonon at the T pointThe nonresonant and resonant 
Raman processes in TMDs will be discussed in Sections 
4 and 5 , respectively. 

The Raman scattering intensity of a Raman-active 
mode in a crystal can be expressed by the Raman tensor 
( R ) in the crystal as da j dXl = \e s • R • e^| 2 , where and 
e s are the polarization vectors of the incoming and scat¬ 
tered photons, respectively. The above equation results 
in the Raman selection rule of R 7^ 0 , where R^ is 
one component of R. Furthermore, \e s -R-ei\ 2 establishes 
the selection rule for the Raman spectrum of the Raman 
mode on and e s . The R value for a Raman mode 
with specific symmetry can be found in some standard 
references E2lLU 


For example, the Raman active modes in bulk MX2 
are denoted as Ai g , E \ g , and 2 E^ g (E^g and E\^). Their 
R values are as follows: 


Alg 


E\g 


E 2g 


fa 0 0 

0 a 0 
0 i 

0 0 0 
0 0c 
^ 0 c 0 

0 d 0 
dO 0 
0 0 


/ 0 0 —c\ 

, 0 0 0 , 

-c 0 0 / 

/d0 0 \ 

, 0 -d 0 • 

y 0 0 0 J 


Under the back-scattering configuration, and e s 
are within the xy plane, which can be set to = 
(cos#, sin#, 0) and e s = (1, 0, 0), where # is the angle be¬ 
tween ei and e s . Usually, ( 1 , 0 , 0 ), ( 0 , 1 , 0 ), and ( 0 , 0 , 1 ) 
of ei and e s are denoted as X, Y, and Z, respectively. 
Thus, the intensity of the A\ g modes (l(Ai g )) in bulk 
MX2 is proportional to a 2 cos 2 #, which monotonically 
decreases from the maximum to the minimum when # 
is varied from 0 {XX) to tt/ 2 {YX). E\ g is doubly de¬ 
generate and there are two Raman tensors: R x and R y . 
l{Ei g )=\e s • R x • e^| 2 + \e s • R y • e^| 2 = 0 , which means 
that it cannot be detected under this polarization config¬ 
uration. Similarly, I(^ 2 p) =d 2 , which is independent of #. 
The calculated Raman spectra of bulk M0S2 in Fig. [ 3 ja) 
show decreasing Ia 1q and constant Ie 29 when # decreases 
from 0 to 7T/2P 

For Raman active A\ g and E g in ENL-MX2, their R 
are listed as follows: 


Alg 


fa 0 0 

0 a 0 

0 i 


E n 


c 0 0 

0 —c d 
\0 d 0 


^0 —c—d\ 

-c 0 0 

\-d 0 0 J 


Ai g in ENL-MX2 has an identical R value to that in 
bulk MX2, indicating similar polarization properties, as 
discussed above. That is, A\ g is absent in the YX con¬ 
figuration and approaches the maximum intensity in the 
XX configuration. However, E g is present in both XX 
and YX configurations. 

R of Raman active A l5 E and E in ONL-MX2 are 
listed as follows: 































7 


E 


/0d0 



E 



In fact, the R value of A ly E , and E are the same as 
those of Ai g , E\ g , and E2 g , respectively. This is because 
of the correlation between D^h and Thus, these 

modes show similar polarization dependence to their cor¬ 
responding modes in bulk MX 2 , as shown in Fig. [3^b) 
for IL-M 0 S 2 . 

Based on the above discussion, we can assign all of the 
observed Raman modes in ONL- and ENL-MX 2 . As an 
example, the Raman spectra of 5L, 6L, and bulk M 0 S 2 in 
the XX and YX configurations are shown in Fig. [3jc). 
Modes that are present in the XX configuration but not 
in the YX configuration correspond to the relative mo¬ 
tions of the atoms perpendicular to the layer plane, and 
are assigned as A 1 (ONL) and A\ g (ENL and bulk). Con¬ 
versely, the modes that exist in both XX and YX cor¬ 
respond to the relative motions of the atoms within the 
xy plane, and are assigned as E (ONL), E g (ENL), and 
E^ g (bulk). The upper right corner of the mode nota¬ 
tion in Fig. [3jc) is used to distinguish modes with the 
same symmetry. Furthermore, a polarized Raman mea¬ 
surement can distinguish two overlapping modes in the 
unpolarized Raman spectrum by their different polariza¬ 
tion dependence. For example, in 5 L-M 0 S 2 , the mode at 
19 cm -1 [E ) can be identified in FJ even though it is 
overlapping with the broad A 1 mode at 17 cm -1 in XX 
. The properties of these modes, including the 2LA(M) 
mode at ^450 cm -1 , will be discussed later. 


IV. NON-RESONANT RAMAN SCATTERING 
IN TMDS 

A. Anomalous frequency trends of E\ g and A± g 
with increasing layer number 

From Table. [I| there are seven Raman active modes 
(Ai g , doubly degenerate ‘lE^g, and doubly degenerate 
Eig) for bulk MX 2 . Under the back-scattering Raman 
configuration, E\ g is absent according to its Raman ten¬ 
sor. Thus, three Raman peaks are expected in the Ra¬ 
man spectrum. Indeed, three peaks at ^33.5 (E^g), 
^383.6 {E^g) and ^408.7 cm -1 (Ai g ) are observed for 
bulk M 0 S 2 , and assigned by their polarization proper¬ 
ties, as shown in Fig. |3|c). The frequencies of the cor¬ 
responding modes in 1L- and ML-MX 2 are expected to 
be dependent on the layer number. Based on the discus¬ 
sion in Section 2, the frequencies of these modes in ONL- 
(E and A x ) and ENL-MX 2 (E g and A\ g ) are different 


to the bulk because of their different symmetries. When 
we discuss the evolution of the peak parameters of the 
Raman modes from the bulk to the monolayer, for the 
high-frequency modes, we simply refer to the correspond¬ 
ing modes as E\ g and A\ g in the bulk, as is commonly 
done in the literature ! 18 * 54 * 

The layer number (N) dependence of the peak position 
and width for E\ g and A\ g of NL-M 0 S 2 is shown in Fig. 
Ufa). Pos(£ 1 2 ^) and Pos(Ai^) show opposite trends with 
(decreasing thickness from the bulk to 1L, as summarized 
in Fig. 0b). However, based on the linear chain model 
(only van der Waals interactions are included), the two 
modes would decrease in frequency from 2L to 1L. 19 This 
unexpected trend of E\ indicates that interactions other 
than van der Waals forces exist which was also re¬ 

vealed by the anomalous Davydov splitting between E\ g 
and E\ u (Section 2). Molina-Sanchez and Wirtz care¬ 
fully examined the self-interaction term of S atoms for 
A\ g (Mo atoms fixed) and both S and Mo atoms for 
E\ g in 1L and bulk MoS 2 .^ They attributed the un¬ 
expected trends of Pos^^) to the long-range Coulomb 
term, which considerably decreases in the FL and bulk 
structures because of the significant increase of the dielec¬ 
tric tensor with NP^ The long-rang Coulomb interaction 
is induced by the effective charges resulting from the rel¬ 
ative displacement between Mo and S atoms. E\ g (Fig. 
|2^a), present in XZ and YZ polarization configurations) 
also shows relative displacement, and is thus expected to 
present an anomalous frequency trend 

Similar anomalous frequency trends with N between 
Ai g and E\ g exist in NL-WS 2 , -WSe 2 , and -MoSe 2 - Fig. 
|4|c) and (d) show the trends for NL-WS 2 P 2 The trends 
for WSe 2 and MoSe 2 are shown in Fig. S2 (ESI). 72 74 
The trends favor the layer number identific ation of MX 2 
through their relative frequency difference.* 18 * 75 * For ex¬ 
ample, the frequency difference Aco(Ai g — E\ g ) of NL- 
M 0 S 2 (N=l-5) can be expressed as Au)(Ai g — E^g) = 
25.8 — 8.4/A, as shown in Fig. 0b). However, fur¬ 
ther N identification for N > 5 must be performed by 
other methods, such as interlayer vibration mode analy¬ 
sis, which will be discussed in Section 4.2. This empirical 
expression accounts for exfoliated NL-M 0 S 2 , but not for 
NL-M 0 S 2 grown by CVD. Au(Ai g — E\^) for CVD-grown 
IL-M 0 S 2 is ~20 cm -1 , which is ^2 cm -1 larger than that 
in exfoliated !L-MoS 2 p^ 


B. Interlayer shear and layer breathing modes 

Because the restoring forces for A\ g jB\ g (LB) and 
E g /E 2 g (C) in ML and bulk MX 2 are weak owing to 
the van der Waals interlayer interactions, the modes are 
expected to have low frequencies. Indeed, the C mode of 
bulk M 0 S 2 is located at 33.5 cm -1 , as shown in Fig. [3^c). 
The LB mode is absent in the Raman spectrum of bulk 
MX 2 because it is Raman inactive. No C or LB modes 
are expected in IL-MX 2 because there is only one X-M- 












FIG. 4: (a) Raman spectra of NL-M 0 S 2 (N=l-8, 10, 14, and 18) and bulk M 0 S 2 . The two grey dashed lines indicate Pos(i? 2 g) 

and Pos(Ai 5 ) in bulk M 0 S 2 . (b) Frequency (uj) of K and Ai g , and the frequency difference between E\ g and A\ g (Ac 0 ) as a 
function of 1/N. For 1 < N < 5, the linear fitting gives Ac o(A — E) = 25.8 — 8.4/A. (c) Polarized Raman spectra of 1L-5L 
and bulk WS 2 , with the frequencies of E\ g and A\ g , frequency difference, l(A\ g )/I(E \ g ), and peak width summarized in (d). 
Reproduced with permission from ref. [72J Copyright 2013, Royal Society of Chemistry. 
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FIG. 5: (a) Stokes and anti-Stokes Raman spectra of NL-M 0 S 2 (N=l-8, 10, 14, and 18) and bulk M 0 S 2 between -60 and 60 

cm -1 . Dashed and dotted lines are guides for the eye for the frequency trends of the C and LB modes, (b) Frequency and 
(c) FWHM of C and LB modes as a function of N. Solid lines are guides for the eye. Reproduced with permission from ref. 
ITqI Copyright 2013, American Physical Society, (d) Diatomic chain model (DCM) and monatomic chain model (MCM) for 
NL-MX 2 , where axx and olxm describe the force constant between X-X layers and X-M layers, respectively, mx and him are 
the mass of the X and M layers, respectively, (e) Symmetry, frequencies, and atomic displacements for the C and LB modes in 
3 L-M 0 S 2 , which are Raman-active (R) and/or infrared-active (IR). 


X layer. However, when IL-MX 2 is vertically stacked to 
form ML-MX 2 , a series of C and LB modes appear in the 
ULF region, whose frequencies will significantly depend 
on the weak interlayer restoring force and layer number 
N. 

The frequencies of LB modes are usually less than 
^100 cm -1 , while those of C modes are smaller (less than 
^50 cm -1 ). The traditional approach to perform very 
low-frequency Raman measurements involves the use of a 
triple spectrometer, which has been applied to ML-M 0 S 2 
or WSe 2 by several groups. 19 20 77 79 A cross-polarized 
backscattering geometry can be used to suppress the 
Rayleigh signal and obtain Raman spectra very close to 


the laser line. 77 Because the spectral throughput is very 
low for the triple spectrometer, the system is difficult to 
apply to various MX 2 compounds and other 2D materi¬ 
als. For example, there have been no reports on the C 
modes in ML graphene by triple spectrometer. Using a 
combination of a single monochromator with BragGrate 
bandpass and notch filters (OptiGrate Corp.), Tan et al. 
detected ULF Raman signals down to ^5 cm -1 with high 
signal throughput. The y success fully observed C and LB 
modes for ML graphene i 17 * 80 * 81 * This ULF technique has 
been applied to multilayer M 0 S 2 up to 18L, where several 
sets of C and LB modes were detected, as shown in Fig. 
[HJa)P 9 Similar results have also been reported for ML 
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Layer number (N) 



Layer number (N) 

FIG. 6: Positions of the (a) C and (b) LB modes of NL-M0S2 
as a function of N. The experimental data are shown as blue 
open circles whose diameter represents the Raman intensity of 
each mode. The green crosses are the calculated results from 
the DCM. The gray dashed lines are another set of the C and 
LB modes based on the MCM. The black solid lines in (a) and 
(b) are fitted by c j(N) = cj(2)^yY^Eos{NoEjN) (N> 2 N 0 ) 
and a; (TV) = cu(2)^/l — cos(TVo7r/TV) (N>2 No), respectively, 
for the phonon branches originating from 2L-M0S2 (No=l), 
6L-M0S2 (No= 3 ), and IOL-M0S2 (No= 5 ) based on the MCM. 
They are simply denoted as and LB^y , respectively. 

Reproduced with permission from ref. [ 19 ] Copyright 2013 , 
American Physical Society. 


WSe 2 up to 7LP Fig. §b) and (c) show the frequency 
and full-width at half maximum (FWHM) of typical C 
and LB modes as a function of TV. The larger FWHM 
of LBM compared with the C mode was interpreted as 
its anharmonic feature with a significant enhancement of 
phonon-phonon scattering with decreasing TV. 79 

To elucidate the origin of these C and LB modes in 
NL-MX 2 , a diatomic chain model (DCM) can be built 
(see Fig. [5jd)). In a DCM, only two force constants are 
required to describe the vibrations under the nearest- 
neighbor interlayer interaction approximation: axx and 
ftXM, where axx is the force constant per unit area 
between the two nearest X planes in two adjacent lay¬ 
ers, and axM is the force constant per unit area be¬ 
tween the nearest X and M planes within a MX 2 layer. 
Their _L (or ||) components are used to describe the vi¬ 


brations along the z axis (or within the xy plane). By 
solving the dynamic eigenequation and fitting the results 
to observed E /A x in 1L-MX 2 and C/LB modes in NL- 
MX 2 (N>2), the strength of intra-/inter-layer coupling 
and the atomic displacements can be obtained. In M 0 S 2 , 

a ss/ a SMo = 2-6% and a ss/ a ^SMo = ^'^°- The atomic dis¬ 
placements of both the C and LB modes in 3L-M0S2 are 
shown in Fig. |5|e), where the corresponding irreducible 
representations are determined by symmetry analysis of 
these displacements. In NL-MX2, there are (TV — 1) LB 
modes and (TV —1) doubly degenerate C modes from sym¬ 
metry analysis (Table. |l|. When N is odd, there are 
E' and e” f or C modes, and A[ and a” 
for LB modes. When N is even, there are y E g and 
E u for C modes, and y A\ g and A 2u for LB modes. 
The peak positions of all of the C and LB modes in NL- 
MX2 can be calculated based on the DCM. Those of 5 L, 
6 L, and bulk M0S2 are summarized in Table. [Il] along 
with the experimental data (Fig. [3jc)) and polarization 
properties. 

Based on the DCM, the relative displacements be¬ 
tween Mo and two S atoms in a S-Mo-S layer are very 
small for all of the C and LB modes in ML-M 0 S 2 . 
For example, the relative displacement in 2 L-M 0 S 2 is 
~ 0 . 6 %, and it decreases with increasing N. Thus, a 
monatomic chain model (MCM) can be considered if the 
entire X-M-X layer is treated as a single ball, as shown 
in Fig. j5jd). Taking the entire layer as a ball with 
mass M = m ,m + 2 mx and interlayer coupling a xx 

for LB modes, and a xx for C modes, the eigenequa¬ 
tion can be analytically solved and the relationship be¬ 
tween the position and N for both the C and LB modes 
can be obtained: 2 N 0 (N) = 2)y/T db cos(TVo7r/TV) 

and uf B , 2 N 0 ( N ) = u lb( 2)-\/l ± cos(V 0 7r/V) (N and N 0 
are integers and TV >2No), where cjc( 2) and cjlb(2) re¬ 
fer to the frequencies of the C and LB modes in 2L- 
MX 2 , respectively. These vibration modes can be as¬ 
signed to a series of branches, as shown in Fig. [6] for 
TVL-M 0 S 2 . Each branch of the C (^C2JV 0 C^O> denoted as 
C^Vq) and LB (u) BB 2 N 0 (N), denoted as LB;^) modes 
includes one high-frequency subbranch (” + ”) and one 
low-frequency subbranch (” — ”). The two subbranches 
cross at TV = 2TV 0 . In the case of the bulk (TV = 
00 ), uj^(bulk) = V2(jUc{2) and uo^ B (bulk) = x/2^lb( 2). 
Thus, the frequencies of the C and LB modes in the bulk 
can be obtained from those detected in 2 L-MX 2 . The 
theoretical results from the DCM and MCM are in good 
agreement with the experimental data for TVL-M 0 S 2 , as 
shown in Fig. [6] Once ujc( 2) and ljlb( 2) are obtained, 
the complete sets of relationships between peak positions 
and TV for both the C and LB modes are completely de¬ 
termined. 

Because cjc( 2) and ujlb{ 2) are linked with the force 
constants between the two nearest S planes in adja¬ 
cent layers (u)c( 2) = (1/ \/2irc)\Jal s /M and ujlb{ 2) = 
(I/V^ttc) y/a-^ 3 /M, respectively), the good agreement 
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TABLE II: Symmetry, polarization (polar.), experimental positions (Exp., in cm -1 ) and theoretical frequencies (Theo., in 
cm -1 ) of the C and LB modes in 5L-, 6L- and bulk M 0 S 2 . The polarization configuration in the bracket is not available under 
the back-scattering condition. 



C 


LB 

5L mode 

E'(R,IR) 

E"(R) 

Ai(R) A 2 (IR) 

polar. 

XX,YX 

(XZ,YZ) 

XX 

Exp. 

19.0, 31.1 

- 

17.5, 46.0 

Theo. 

19.1, 31.0 

- 

17.5, 46.1 

6 L mode 

E 9 (R) 

E„(IR) 

Ais(R) A 2 „(IR) 

polar. XX,YX,(XZ,YZ) 

- 

XX 

Exp. 

22.7, 31.2 

- 

14.6, 40.3 

Theo. 

23.0, 31.5 

- 

14.6, 40.2 

Bulk mode 

El*(R) 


Bl,(S) 

polar. 

XX,YX 


- 

Exp. 

33.5 


- 

Theo. 

33.5 


- 


between the MCM and the experimental data for NL- 
M 0 S 2 is indicative of a constant interlayer interaction 
from 2L to ML to bulk M 0 S 2 I 19 Other TMDs show sim¬ 
ilar structures to M 0 S 2 , which indicates that the above 
equation can also be applied to these materials by replac¬ 
ing the 2) and cjc(2 ) values with the corresponding 
values of the other 2L-TMDs. Furthermore, this MCM 
can predict the dependencies of the C and LB modes on 
the layer thickness in other LMs, such as hexagonal boron 
nitride (hBN)P^It should also be noted that the C and 
LB modes are directly related to the interlayer coupling, 
which has the potential to probe the interface coupling 
in folded (or twisted) TMDs and heterostructures built 
from TMDs. 


C. Phonon dispersion trends away from the T point 



Because of the Heisenberg uncertainty principle, the 
fundamental q ~ 0 Raman selection rule is relaxed for 
a finite size domain, allowing participation of phonons 
away from the Brillouin-zone center (T). The phonon 
uncertainty is A q ^ 1/d, where d is the size of the crys¬ 
tal domain. The Raman spectrum of IL-MX2 with small 
domain size can provide sufficient information about the 
phonon dispersion near the T point. This usually gives 
a downshift and asymmetric broadening of the Raman 
peak. The corresponding Raman mode will usually red- 
shift in frequency.^ 83 The ion-implantation technique 
has been widely used to obtain nanocrystalline graphene, 
ML graphene, and graphite! 84 * 85 * Recently, small crystal 
size IL-M0S2 and IL-WS2 were obtained by increasing 
the dosage of ion implantation. 86 As shown in Fig. [^c) 
andj^d), with increasing ion-implantation dosage, soft¬ 
ening of E and hardening of A x in IL-M0S2 occurs, 
while both the E and A x modes are softened in 1 L- 


FIG. 7: Phonon dispersions of ZO 2 (A 1 at T), LO 2 (E at 
T), and TO 2 (E at T) along T-M for (a) IL-M 0 S 2 and (b) 
IL-WS 2 . Experimental (crosses) and calculated (solid lines) 
results of (c) IL-M 0 S 2 and (d) IL-WS 2 for different ion im¬ 
plantation doses. The dashed lines are the contribution of the 
subpeaks from the LO 2 and TO 2 branches. The spectra are 
normalized to A 1 . 


WS 2 . Correspondingly, ion implantation also results in 
the reduction of the PL intensity from the A exciton. PL 
is almost quenched when the ion-implantation dosage is 
larger than 5xl0 13 Ar + /cm 2 , which indicates that 1L- 
M 0 S 2 is highly disordered. 

To understand the above results, the phenomenolog¬ 
ical phonon confinement model originally proposed by 
Richter, Wang, and Ley (WRL model/ 87 * has been ap¬ 
plied to nanocrystalline IL-M 0 S 2 and IL-WS 2 induced 
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by ion-implantation. This model explains the A[ and 
E frequency shifts of IL-M 0 S 2 with increasing ion- 
implantation dosage, as shown in Fig. [7^c) by vertical 
dash-dotted lines. The small hardening of A 1 in IL-M 0 S 2 
results in the ZO 2 phonon branch slowly increasing in fre¬ 
quency when q is away from the T point, which agrees 
with the trend of ZO 2 dispersion in Fig. [7ja). Simi¬ 
larly, the larger softening of E in IL-M 0 S 2 suggests that 
LO 2 /TO 2 branches rapidly decrease with q away from 
the T point. However, for nanocrystalline IL-WS 2 , both 
A x and E are softened with increasing ion-implantation 
dosage (Fig. [T^d)), which can be well-reproduced by the 
WRL model. When the ion dosage increases, the trend 
of phonon dispersion away from the Y point is steeper, 
the peak frequency shift is larger, and the peak profile is 
more asymmetric and broader. Although LO 2 of IL-WS 2 
rapidly increases near the M point (Fig. (7^b)), only non¬ 
zero q close to the Y point can participate in the Raman 
process, so the corresponding blueshift of E cannot be 
observed. 


D. Raman characterization of TMD alloys 



FIG. 8: (a) Raman spectra of Moi-ccWa^ monolayers with 
different W composition x. The three solid lines are guides 
for the eye to show Pos(E ) and Pos(A 1 ) peaks with W com¬ 
position x. (b) Composition-dependent Raman frequencies 
of each mode in Moi-^Wa^ monolayers. Reproduced with 
permission from ref. [88] Copyright 2014, Royal Society of 
Chemistry. 


Alloying materials with different bandgaps have been 
widely used in bandgap engineering of bulk semiconduc¬ 
tors. The realization of bandgap tuning in atomically 
thin 2D alloys with various choices of layered TMDs could 
allow a wide range of bandgap tuning for application 
of 2D materials in nanoelectronics and optoelectronics. 
Monolayer TMD alloys have been proposed by theor etical 
calculation^® and demonstrated in experiments 9192 . 
A series of monolayer TMD alloys (Moi_ a ,W ;E S 2 , x = 
0 to 1) have been cleaved from their bulk crystals. 91 A 
monolayer alloy film of MoS 2 (i- cc )Se 2 CC of arbitrary com¬ 


position can be grown by CVD through control of the 
S/Se ratio of the organic precursors used in the growth 
process.® 

Structure characterization of 2D alloys, such as com¬ 
position and atom mixing, is of fundamental impor¬ 
tance for their potential applications. Raman spec¬ 
troscopy is a powerful tool to characterize structures of 
2D alloys by the frequency shift for alloy composition 
and by peak broadening for alloy degree.® Compared 
with FL and bulk Moi_ a ,W ;E S 2 , Moi-^W^^ monolayers 
show a smaller frequency difference between A 1 and E 
and similar polarization dependence of Raman modes.® 
Composition-dependent unpolarized Raman spectra of 
Moi_ a ,W a ,S 2 (x = 0-1) monolayers are shown in Fig. 
[ija). The assignments of A 1 and E are based on the 
polarized Raman spectra? 9 Similar to Moi_ a ,W cc S 2 bulk 
alloys,® A x and E show one-mode and two-mode behav¬ 
ior, respectively. As shown in Fig. Ib), with increasing 
x from 0 to 1, the A ± mode continuously upshifts, while 
the MoS 2 -like E mode shifts to lower frequency. WS 2 - 
like E appears at x = 0.2 and slowly moves towards 
higher frequency. As the W composition increases, the 
relative intensity of WS 2 -like E increases, while that 
of MoS 2 -like E decreases. The Raman peak at ^360 
cm -1 present from x = 0.13-0.77 is assigned to an al¬ 
loy disorder-related peak, which is not sensitive to the 
alloy composition. The peak frequency at ^346 cm -1 
is attributed to 2LA(M) of WS 2 , and overall shows a 
decreasing trend with decreasing x (Fig. [8^b)). Based 
on the modified random-element-isodisplacement model, 
Chen et al. reproduced the dependence of A 1 , MoS 2 -like 
E , and WS 2 -like E on x, which can be used to quantify 
Mo/W compositions.®- 95 This systematic Raman scat¬ 
tering investigation of Moi_ cc W ;e S 2 monolayers is very 
helpful for the spectroscopic determination of the com¬ 
position and alloy degree of various 2D crystal alloys. 


V. RESONANT RAMAN SCATTERING IN 
TMDS 

A. Phonons of bulk MX 2 at the Y and M points 

First, we will discuss the nonresonant and reso nant 
Raman spectra of bulk M 0 S 2 , as shown in Fig. [ 9 ]® The 
nonresonant and resonant Raman spectra of bulk WS 2 
and WSe 2 are shown in Fig. S4 (ESI). The nonresonant 
spectrum is usually dominated by four basic vibration 
modes lying at the Y point EEEMeEH E% g (32 cm -1 , C), 
Ei g (286 cm -1 ), E\ g (384 cm -1 ), and A\ g (409 cm -1 ). 
Because E\ g (286 cm -1 ) is forbidden under the backscat- 
tering geometry according to its Raman tensor, it does 
not appear in the spectrum excited by 2.33 eV in Fig. 
Ufa). A very weak peak at ^450 cm 1 (2LA(M)) and 
several weaker peaks above 520 cm -1 are also observed 
in the nonresonant Raman spectrum in Fig. |9ja). 

Because 1.96 eV is close to the energy of the A exci- 
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FIG. 9: (a) Nonresonant and (b) resonant Raman spectra 

of bulk M 0 S 2 at room temperature. The peak position and 
corresponding assignment of each mode are shown, where T, 
M, and T -A in brackets after each mode indicate the involved 
phonons that come from T, M points, and along the T-A 
direction, respectively. The asterisk (*) labels the Raman 
peak of Si —520 cm -1 . Reproduced with permission from ref. 
l96l Copyright 2014, AIP Publishing LLC. 


ton in bulk M 0 S 2 (~1.88 eV at room temperature) P it 
is usually used to obtain the resonant Raman spectrum 
of bulk M 0 S 2 , as shown in Fig. [9jb)P^ Under resonant 
excitation, intense modes are found atP^ 7 ^~178 cm -1 
(A lg (M)-LA(M)), -420 cm” 1 (B* g +E* u (r-A)), -465 
cm -1 (2LA(M) (—450 cm -1 ) and A 2n (r) (—465 cm -1 )), 
-572 cm -1 (2E lg (T)), -600 cm -1 [E\ g (M)+ZA(M)), 
and —643 cm -1 (Ai g (M)+LA(M)). The mode at —420 
cm -1 is highly dispersive and explained by Sekine et al. 
as a second-order Raman process, where a polariton cre¬ 
ated by the incident laser in a photon-like branch will 
first be scattered by emitting a longitudinal dispersive 
quasiacoustic (QA) phonon along T-A, and then scat¬ 
tered by emitting a E\ u phonon (which is dispersion- 
less along T-A) into the exciton-like final stated Here, 
T-A is perpendicular to the 2D Brillouin plane. The 
QA phonon at the T point is B\ of M 0 S 2 . The en¬ 
ergy difference of its Stokes and anti-Stokes shifts un¬ 
der 217 K is —4.7 cm _1 P Pos(H|^) can be deduced 
from the dispersive B^gEE^iY-A) based on the exciton- 
polariton dispersion.^ The most intense and asymmet¬ 
ric resonant mode at —465 cm -1 was first assigned by 
Stacy et al. as 2LA(M)pSI an d then suggested by Frey 
et al. to be composed of a double-mode feature (second- 
order 2LA(M) and first-order A 2u (r)), because of its 
splitting in nanoparticlesThe mode at —643 cm -1 is 
assigned as Ai g (M)+LA(M), which can be used to es¬ 
timate the energies of Ai g (M) and LA( M) when com- 
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FIG. 10: Resonant Raman spectra of bulk M 0 S 2 at room 
temperature (top) and at 4.2 K (bottom). The black arrows 
indicate the five phonon replicas of Ei g (M) and A 2 u(M). The 
mode at 782 cm -1 is the first phonon replica of the mode at 
603 cm -1 , as indicated by the blue line with double arrows. 
Reproduced with permission from ref. 1961 Copyright 2014, 
AIP Publishing LLC. 


bined with the mode at 178 cm -1 (Ai g (M)-LA(M)). 
Thus, A lg (M) and LA{ M) are at —410.5 cm -1 and 
—232.5 cm -1 , respectively.^ The frequency of Ai g (M) 
(—410.5 cm -1 ) is close to that of Ai q(r) (—409 cm -1 ), 
reflecting its dispersionless character ! 53 l 54 l The mode at 
—600 cm -1 is £ 1 2 ^(M)+LA(M), which leads to E\ g { M) 
at —367 cm -1 . Similar to Ai g (M)-LA(M) (—178 cm -1 ), 
the mode E 2 g (M)-LA(M) is expected at —134 cm -1 . 
However, the obtained LA(M) (—232.5 cm -1 ) is larger 
than the value determined from 2LA(M) in the nonres¬ 
onant spectrum (—226.3 cm -1 ). Recently, the double 
structure of the mode at —465 cm -1 was reassigned by 
Golasa et al. as 2LA(M) and the first phonon replica 
(Ei g (M)+XA(M)) of Ei g (M) inspired by the close anal¬ 
ogy in the profile to the mode at —643 cm -1 . The mode 
at —643 cm -1 is accordingly assigned as Ai g (M)+LA(M) 
and the second phonon replica ( Ei g (M)+2XA(M )) of 
E'ip(M), where XA is related to the transverse acoustic 
(TA) and/or out-of-plane acoustic (ZA) modeP^ 

The Raman spectrum at low temperature can be used 
to identify the two-phonon difference combination mode. 
Fig. [lO] shows the Raman spectra of bulk M 0 S 2 at 
room temperature (upper panel) and 4 K (lower panel). 
The intensity of a two-phonon mode at — UJ 2 will be 
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highly dependent on the temperatures for creation of 
phonon and annihilation of oui:n(uJi,T)[n(uj 2 ,T) + 1], 
where n(u;,T) = (e huj / kT — l) -1 is the phonon occupa¬ 
tion probability.^" The intensity of the mode at —178 
cm -1 at 4 K is predicted to be 37 times weaker than 
the intensity at room temperaturei 9 ^ Thus, this peak 
disappears at low temperature, as shown in Fig. [10] 
That is why the mode at —178 cm -1 is assigned as 
Ai g (M)-LA(M) in Fig. j9^b). Moreover, Golasa et al. 
observed five phonon replicas of E\ g { M) (—466, —645, 
-825, -1005, and -1185 cm- 1 ) and A 2u { M) (-574, 
—751, —930, —1108, and —1290 cm -1 ) at 4.2 K, as shown 
in Fig. [To] (lower panel). In addition, the mode at —782 
cm -1 (E 2 g(M)+LA(M)+XA(M)) can be considered as 
the first phonon replica (E 2 g (M)-\-LA(M)) of the mode 
at —603 cm -1 , as shown by blue line with double ar¬ 
rows in Fig. [lO] (lower panel). This leads to a series of 
phonon replicas at —962 (E2 g (M)+LA(M)+2XA(M)), 
1142 cm- 1 ((£ ,1 ^(M)+LA(M)+3XA(M))), and so on. 
The resulting frequency of LA(M) is —230.5 cm -1 , and 
XA(M) is —180 cm -1 . LA(M) approaches the value de¬ 
termined from the nonresonant spectrum, but the con¬ 
tribution of the involved XA phonon to the M point re¬ 
quires further theoretical analysis.^ From the phonon 
dispersion in ref. [54] XA probably corresponds to the 
ZA branch (See Fig. |2jc)). 

Other bulk MX 2 compounds are expected to exhibit 
similar resonant behavior to bulk M 0 S 2 . However, the 
resonant excitation energy will be different because of 
the different band structures. For example, Fan et al. 
found that the Raman resonance in bulk WS 2 covers a 
large energy range from 1.96 (633 nm) to 2.41 (514 nm) 
eVjlQU which is attributed to the broad features of its B 
exciton (Fig. S4, ESI) 


B. Phonons of monolayer MX 2 at the T and M 
points 


Raman investigation of IL-MX 2 is critical to under¬ 
stand the phonon spectrum of FL-MX 2 . Fig. 11 shows 
the resonant and nonresonant Raman spectra of 1L- 
M 0 S 2 , lL-WSe 2 , and IL-WS 2 . The insets show the de¬ 
tailed fitting results for the resonant peaks near A 1 and 


E . The peak notations in bulk MX 2 have been changed 
in 1L-MX 2 because of the different point group (Dsh) of 


1L-MX 2 . 

In the case of the Raman spectrum of IL-M 0 S 2 in 
Fig. [llja), the modes at —178, —450, and —628 cm -1 
are assigned as H. / 1 (M)-LA(M), 2LA(M), and (M)- 
LA(M), respectively, which is similar to the bulk. 9 ^- 00 * 1 ^ 3 ! 
Similarly, the modes at —150 and —600 cm -1 are as- 
signed as E (M) l ° 2 -LA(M) and E (M) L ° 2 +LA(M) PTH 
Here, E (M) LO2 is the high-frequency branch of the two 
split branches of E at the M point (E'(M) L ° 2 and 
E (M) t ° 2 ), as shown in Fig. [^b). The mode at —600 
cm -1 is E (M) l ° 2 +LA(M), with its second and third 





Raman Shift (cm' 1 ) 


FIG. 11: Resonant (black solid lines) and nonresonant (gray 
dashed lines) Raman spectra of (a) IL-M 0 S 2 (b) lL-WSe 2 (c) 
IL-WS 2 at room temperature. The insets show the subpeak 
fitting results. The frequencies of each mode are shown. The 
asterisk (*) indicates the Raman peak of Si at —520 cm -1 . 


phonon replicas found at —638 (E (M) TOl +2ZA(M)) 
and 820 (E (M) TOl +3ZA(M)) cm -1 , respectively. The 
modes at —756, —782, and —820 cm -1 are 2E' (M) l ° 2 , 
E (M) l ° 2 +A 1 (M), and 2A 1 (M), respectively, which can 
be used to estim ate the spin-orbit splitting in IL-M 0 S 2 
(see Section 5.4). 104 l It should be noted that the mode 
at —420 cm -1 (B 2 g -\-Ef u (T-A)) in bulk M 0 S 2 disappears 
because of the absence of the interlayer breathing mode 
in IL-M 0 S 2 . 

For the Raman spectrum of lL-WSe 2 in Fig. 
[n];b), E' (M) l ° 2 -LA(M), y(M)-LA(M), 2LA(M), 

E(M) i02 +LA(M), and A' 1 (M)+LA(M) are found at 
—116, —136, —261, —374, and —398 cm -1 , respectively. 
The modes at —96 and —157 cm -1 are F/(M) t ° 2 -LA(M) 
and E(M) l ° 2 -TA(M), where E (M) T ° 2 +LA(M) (ex¬ 
pected at —357 cm -1 ) and E\ M) l ° 2 +TA(M) (expected 
at —332 cm -1 ) are not observed. The combination 
modes of LA(M), TA(M), and ZA(M) are observed at 
-219 (LA(M)-I-TA(M)), -239 (LA(M)+ZA(M)), -263 
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(3TA(M)), and —394 cm -1 (3LA(M)). The modes at 
—634 and —656 cm -1 are E\ M) l ° 2 +3LA(M) and 
A 1 (M)+3LA(M), respectively. E (M) l ° 2 +2LA(M) and 
A / 1 (M)+2LA(M) are expected at —504 and —528 cm -1 , 
but they are easily distinguished because of the broad 
peak range from 460 to 540 cm -1 in Fig. 


11 


b). 


For the Raman spectrum of IL-WS 2 in Fig. |ll[c), 

#'(m) t ° 2 -la(m), f/(m) l ° 2 -la(m), a / 1 (m)-la(m), 

2LA(M), and Ai(M)+LA(M) are at -176, -193, -231, 
—353, and —583 cm -1 , respectively.^ In addition, the 
modes at —204 and —214 cm -1 can be assigned as 
E"( M) TOi -ZA(M) and E” (M5 TGi -TA(M). 4LA(M) is 
also observed at —702 cm -1 p2I 

All of the tentative assignments of the main peaks ob¬ 
served in IL-M 0 S 2 , lL-WSe 2 , and IL-WS 2 are summa¬ 
rized in Table. uni 


C. Layer-dependent resonant Raman scattering in 
few layer M 0 S 2 



FIG. 12: Resonant Raman spectra of NL-M 0 S 2 (N = 1, 2, 
4 and 7) and bulk M 0 S 2 excited by a 633 nm laser. Different 
modes are indicated with arrows. The symmetry assignments 
are also given for each mode. A noticeable feature is the 
appearance of completely new modes and their evolution with 
layer number. Reproduced with permission from ref. 11051 
Copyright 2013, John Wiley & Sons, Inc. 


different notations of the modes in ENL, ONL, and bulk 
M 0 S 2 are ignored for the convenience of comparison. For 
the mode assignments in bulk M 0 S 2 , refer to Section 5.1, 
except for the observation of Bi u (T) at —409 cm -1 and 
£’i p (M)+LA(M) at —526 cm -1 . In the low-frequency 
spectral region (160-250 cm -1 ), the lineshape asymmetry 
of the Ai p (M)-LA(M) mode at —179 cm -1 is more dis¬ 
tinct for N = 2, 4, and 7 because the A\ g branch is nearly 
flat and the LA branch is concave at the M point (see Fig. 
|2|b) and (c)). The additional peak near 230 cm -1 for N 
= 2, 4, and 7 can be assigned as the first-order LA(M). 
In the spectral region (360-430 cm -1 ), the mode at —377 
cm -1 that is broad in IL-M 0 S 2 becomes sharper in 2L-, 
4L- ,and 7 L-M 0 S 2 . It is assigned as E\ u with a finite 
wavevector near the T point along T-A (c axis) which 
is in contrast to the previous assignment (Fq n (r))! 98 * 99 * 
Actually, E? u { T) is observed on the higher frequency side 
of E , 2 ^(r)Plts relative integrated intensity with respect 
to E^giV) decreases with increasing N, as shown in Fig. 
12 The prominent peak at —420 cm -1 is assigned as 


B^gEEfn (T-A), which can be used to approximately pre¬ 
dict the frequency of B^giT- A) as —43 cm -1 . Bi u (T) at 
—403 cm -1 is the Davydov couple of Ai g (T) (see Fig. 
Ha)), and its relative intensity with respect to A\ g in¬ 
creases with increasing N. The frequencies of 2LA(M) 
and A 2u (T) increase with increasing N, but the relative 
intensity of 2LA(M) to A 2u ( T) decreases with increasing 
N, which is similar to the intensity ratio dependence on 
the size of nanoparticlesThe spectral region between 
430 and 490 cm -1 is assigned to the second-order Raman 
spectra, which also shows layer dependence with high in¬ 
tensity. The modes ,Ei^(M)+LA(M), E^(M)+LA(M), 
and Ai^(M)+LA(M) have an asymmetrical tail towards 
the low-frequency region, which is because of the nearly 
flat dispersion for E\ g and A\ g and concave down LA 
along T-M. There are several new peaks at —554 (Q), 
—588 (P), and —628 (S) cm -1 in 2L-, 4L-, and 7 L-M 0 S 2 , 
whose peak positions do not vary with N, although their 
relative intensities with respect to the neighboring modes 
clearly depend on N.^H A quantitative understanding of 
these new features would require theoretical calculations 
of phonon dispersion and two-phonon density of states as 
a function of N. 


The number of atoms in the unit cell of NL-LMs in¬ 
creases with increasing N. More vibrational modes are 
expected to be involved in the Raman scattering as N 
increases. Because of the existence of weak interlayer 
coupling in 2D LMs, FL-LMs show a distinctive band 
structure, resulting in the dependency of the correspond¬ 
ing interband transition energy on N. Thus, resonant Ra¬ 
man scattering in 2D LMs is expected to be layer-number 
dependent. 

Fig. [T2| shows the resonant Raman spectra of NL-M 0 S 2 
(N = 1, 2, 4, and 7) and bulk MoS 2 -^^The strong Raman 
mode marked with an asterisk (*) is from the Si 02 /Si 
substrate and is completely absent in bulk M 0 S 2 . The 


D. Spin-orbit splitting in M 0 S 2 revealed by triply 
resonant Raman scattering 

The sizeable spin-orbit (SO) splitting in IL-M 0 S 2 is 
promising for its potential application in spintronics and 
quantum information processing, but also makes it a 
challenge to develop an experimental method to detect 
and explore the intrinsic SOC. Previous attempts have 
mainly focused on the energy difference between the A 
and B exciton peaks in the correspo nding one-photon 
(linearand two-photon (non-linear) 40 * 41 * ^ * PL spectra. 
However, the identifications of the A and B exciton peaks 
are still uncertain because of the involvement of an addi- 
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TABLE III: Peak positions and corresponding assignments of each Raman modes observed in IL-M 0 S 2 , lL-WSe 2 and IL-WS 2 . 


1 L-M 0 S 2 

Peak (cm -1 ) 

Assignments 

lL-WSe 2 

Peak (cm -1 ) 

Assignments 

IL-WS 2 

Peak (cm -1 ) 

Assignments 

150 

F , (M) z '° 2 -LA(M) a 

96 

f'(M) t ° 2 -LA(M) 

176 

f'(M) t ° 2 -LA(M) 

178 

Ai(M)-LA(M) a 

116 

E (M) L ° 2 -LA(M) b 


or LA(M) d 

376 

e\M) l ° 2 

136 

A'i(M)-LA(M) 

193 

E (M) b ° 2 -LA(M) b 

385 

F(r) a 

157 

f'(M) z '° 2 -TA(M) c 

204 

F"(M) TOl -ZA(M) e 

404 

4( ry 

219 

LA(M)+TA(M) C 

214 

F"(M) TOl -TA(M) e 

414 

4(m) 

239 

LA(M)+ZA(M) C 

231 

Ai(M)-LA(M) d 

450 

2LA(M) a 

248 

e'( r) b 

298 

E"(T) 

460 

F"(M) TOl +ZA(M) 

250 

A[( r) b 

325 

e"( M) TOie 

566 

2 s"(r) a 

258 

A2(M) C 

347 

E { M) T ° 2e 

600 

S(M) z '° 2 +LA(M) a 

261 

2LA(M) b 

353 

2LA(M) d 

628 

A 1 (M)+LA(M) a 

263 

3TA(M)° 

357 

F(r) d 

638 

E"{ M) TOi +2ZA(M) 

374 

£'(M) i02 +LA(M) b 

419 

a[( r) d 

756 

2S'(M) i02a 

394 

3LA(M) b 

583 

Ai(M)+LA(M) d 

782 

E (M) l ° 2 +A 1 (M) a 

398 

A' 1 (M)+LA(M) 

702 

4LA(M) d 

820 

e" (M) tc>1 +3ZA(M) 

534 

2A' 1 (M) 




2Ai(M) a 

634 

£'(M) iC,2 +3LA(M) 





656 

Ai (M)+3LA(M) 



a refH03l 

b ref [72, c based on the phonon dispersion in 

refl63l d refl65l e based 

on the phonon dispersion in ref!54l 




360 400 760 800 840 

Raman Shift (cm' 1 ) 


’ ■ if 1 

488 nm 


1L 

Av. 

2L 

i_A 

4L 

AA 

bulk 

360 400 760 8(H) 840 


Raman Shift (cm' 1 ) 



FIG. 13: Raman spectra of 1L, 2L, 4L, and bulk M 0 S 2 
excited by (a) 325, (b) 457, and (c) 488 nm laser lines. All 
of the spectra are normalized by I(Ai g ). The dashed lines in 
(a) are the fitting curves. The vertical dashed lines represent 
the peak positions of the Raman modes in bulk M 0 S 2 . (d) 
Schematic diagram of the TRRS process in M 0 S 2 , where 
and V 2 represent the valence band splitting because of the 
spin-orbit interaction, and ci represents a conduction band. 
Reproduced with permission from ref. 11041 Copyright 2013, 
American Physical Society. 


tional bound exciton peak that is sensitive to the ambient 
environment.^ 4 106 The electron (hole) in the CB (VB) 


band resonantly excited by incident laser can be scattered 
into a real electronic state by emitting a phonon, result¬ 
ing in a double-resonant (DR ) or triply resonant (TR) 
Raman scattering processThese resonant Raman pro¬ 
cesses commonly exist in graphene. The resonant Raman 
spectrum is thus very sensitive to the matching between 
the VB splitting energy and the energy of the phonons 
involved. 

Fig. [T3^a), (b), and (c) show the Raman spectra of 1L, 
2L, 4L, and bulk M 0 S 2 excited by 325, 457, and 488 nm 
laser lines, respectively. 10 ^ The three high-order modes in 
the range 750-840 cm -1 excited by 325 nm laser are dra¬ 
matically enhanced compared with the other laser lines, 
and are assigned as the second-order combination Ra¬ 
man modes of E\ g and A\ g . These three Raman modes 
are also present in IL-M 0 S 2 at ~756, ^782, and ^820 
cm -1 (Fig. [lT|^a)). The intensity enhancement has been 
attributed to the electron-two-phonon involved TR Ra¬ 
man scattering (TRRS) process. In the TRRS process 
(Fig. [Tap)), (1) an electron is excited from the v 2 band 
to the Ci band (cq band in ref. 104 at the bandgap) by 
absorbing a photon at the K point of the Brillouin zone, 
(2) the hole in the v 2 band is scattered by a phonon with 
momentum q to the v\ band by an interband transition 
as a consequence of the deformation potential interac¬ 
tion, (3) another phonon with momentum -q scatters the 
hole to the top of the v\ band by intraband transition to 
form an exciton with the e lectr on in the band because 
of the Frohlich interaction^^ and (4) the electron-hole 
pair recombines at the top of the v\ band and emitting 
a photon. For energy and momentum conservation in 
the TRRS process, Ei aser — E exc n on T ‘ZE^i , Ei aser — 

2gr 
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Eexciton + 2 E A lg , and Elaser — -^exciton + E E\ g + E A lg 

for the 2 E\ g , 2Ai g , and E\ g + A\ g modes, respectively. 
The flat phonon dispersion favors energy and momentum 
conservation, where the correct q value can be selected in 
a large range along T-M for the correct phonon energy 
(AEso/ 2), which is responsible for observation of the 
strong TRRS. In experiments, the three resonant high- 
order peaks in IL-M0S2 excited by 325 nm laser have en¬ 
ergies between 93 ( 2 E\ g , which is actually 2 E ) and 102 
meV ( 2 Ai g , which is actually 2A 1 ), giving a SO splitting 
of ^100 meV for IL-M0S2 A^It should be noted that the 
splitting in 2L-M0S2 arises from the combination of SO 
coupling and interlayer couplingP Much weaker intensi¬ 
ties of overtone and combination modes are observed in 
bulk M0S2 (Fig- p^a)). Bulk M0S2 is predicted to have 
much larger VB splitting, which results in a decrease in 
the probability of TRRS because of the loss of energy 
match between VB splitting and the energy of the two 
phonons involved. 


VI. RAMAN SCATTERING OF TMDS 
MODIFIED BY EXTERNAL PERTURBATION 


External perturbations can modify the electric, 
phonon, thermal^ and mechanical properties of 1 L-, 2 L-, 
and ML-TMDs, Chilli which facilitates their application 
in tunable photonic devices, solar cells, flexible electron¬ 
ics, thermoelectric energy conversion, field effect transis¬ 
tors, and catalysts for hydrogen evolution 11 ^. For ex¬ 
ample, bandgap engineering throug h str ain has been re¬ 
ported for both 1L- and 2L-M0S2 p dSHH D Chang et al cal¬ 
culated the electronic properties of M0S2, MoSe2, WS2, 
and WSe2 monolayers, and found that their bandgaps are 
more sensitive to biaxial strain than uniaxial strain. 112 
A large value of the Seebeck coefficient has been re¬ 
ported for IL-M0S2 in thermoelectric energy conversion 
(- 4 xl 0 2 and -lxlO 5 /iVK -1 depending on the strength 
of the external electric field). 114 Electrical control of 
neutral and charged excito ns in several lL-TMDs has 
been recently reported) 117 118 Mitioglu et al. observed PL 
peaks associa ted with charged excitons of IL-WS2 at low 
temperaturesP 19 The pressure-induced phase transition 
is revealed in M0S2P- 20 TMDs under external perturba¬ 
tion has also attracted considerable attention. 

Raman spectroscopy is an ideal technique to probe the 
influence of external environmental conditions on mate¬ 
rial systems by spectral-feature analysis with high reso¬ 
lution. This makes Raman spectroscopy useful in funda¬ 
mental research and device characterization. Indeed, the 
phonon spectra of 1L- and ML-TMDs are significantly 
affected by external perturbation, and their response can 
be probed by Raman spectroscopy. In the following, we 
will discuss how external perturbations such as strain, 
pressure, electric field, charge transfer, temperature, and 
substrate affect the Raman spectra of TMDs. 


A. Effect of uniaxial and biaxial strain 


Biaxial and uniaxial strains can be applied to a sam¬ 
ple. Sahin et al recently calculated the phonon disper¬ 
sion of lL-WSe2 after 1% biaxial compressive and tensile 
strain was applied, and found that such biaxial defor¬ 
mation only results in collective softening/hardening of 
the vibrational modes and does n ot lif t the degeneracy of 
the E and A 1 modes of lL-WSe2 Chang et al. calcu¬ 
lated the phonon dispersions of M0S2, WS2, WSe2, and 
MoSe2 monolayers, and found that their phonon frequen¬ 
cies shifted upward/downward under compressive/tensile 
biaxial strain.^ Fig. 14 'a) and (b) show the phonon dis¬ 


persion of IL-M0S2 and IL-WS2 under no strain, and 
compressive (biaxial), tensile (biaxial), armchair (uniax¬ 
ial), and zigzag (uniaxial) strain. E and E are degen¬ 
erate under both compressive and tensile strain, which 
indicates that these modes remain degenerate as long as 
the hexagonal symmetry of MX2 is retained. However, 
these two modes split when uniaxial strain is applied, in¬ 
dicating that the isotropic symmetry in the xy plane has 
been broken. 


A three-point bending apparatus (Fig. [l 4 [c)) is usu¬ 
ally used to apply uniaxial strain to small MX2 flakes 
Zhu et al observed the splitting of E /E g with increasing 
strength of uniaxial strain in both 1 L- and 2L-MoS2p^ 
as shown in Fig. 


14 


d) and (e). A 1 /A\ g describes the 

vibrations along the z axis, whereas E /E g describes the 
vibrations in the xy base plane. In the absence of strain, 
E /E g are 2D degenerate modes. Applied uniaxial strain 
will remove the degeneracy because of symmetry break¬ 
ing within the xy plane, leading to mode splitting for 
E and E g . The lattice symmetry of 1L- and 2L-M0S2 
changes to C s and C2/1 under uniaxial strain, respectively. 
The corresponding mod es a re characterized by three A 
(1L) and three A g (2L)P^ For IL-M0S2, with increas¬ 
ing applied strain, the frequency shift is small for the 
A ( 3 ) mode, indicating that the strain has little effect on 
vibrations perpendicular to the xy plane (Fig. Eld)). 
The larger frequency shift of A ( 1 ) with increasing strain 
indicates that it vibrates along the direction of the ap¬ 
plied strain. The intensities of the split A ( 1) and A (2) 
modes orthogonally respond to the angle between the po¬ 
larization of the scattered Raman signal and the strain 
axis, which can be used to identify crystallographic ori¬ 
entations of lL-MoS2f 122 Similar results have also been 
observed in uniaxially strained 2L-M0S2 (Fig. [T^e)). 
Moreover, the dispersion of low-energy flexural phonons 
in strained lL-MoSe2 has been found to cross over from 
quadratic to linear.^ Based on the strain dependence of 
A (l)and A (2), Conley et al determined parameters to 
characterize the anharmonicity of the molecular poten¬ 
tials: the Grimeisen parameter (7) and the shear defor¬ 
mation potential (J 3 ). 7^ = 1.1 ± 0.2 and f 3 E > = 0.78 =b 
0.1 were obtained. 11 ^ In addition, strain naturally exists 
in state-of-the-art TMD nanotubes, which is similar to 
MX2 with uniaxial strain. 123 
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FIG. 14: Phonon dispersion of (a) IL-M 0 S 2 and (b) IL-WS 2 . For compressive (tensile) biaxial strain, the phonon frequencies 
of E and E shift upward (downward). For uniaxial strains, the degeneracies of E and E at T lift because of the introduced 
anisotropy of bond strengths. Reproduced with permission from ref. 1661 Copyright 2013, American Physical Society, (c) 
Three-point bending apparatus, (d) E splits into two modes as tensile strain is applied to IL-M 0 S 2 . (e) Same as (d) but for 
2 L-M 0 S 2 . Reproduced with permission from ref. 11161 Copyright 2013, American Physical Society. 


B. Temperature dependence 


TABLE IV: Comparison of the temperature coefficients (xt) 
for each mode in 1L-, FL- and bulk TMDs. 


Material 

Type 

El/E' 

(cm 1 /K) 

Alg/Al 

(cm 1 /K) 

IL-M 0 S 2 

silicon-supported 

-0.013 a 

-0.016 a 

IL-M 0 S 2 

suspended 

-0.011 b 

-0.013 b 

IL-M 0 S 2 

sapphire-supported 

-0.017 b 

-0.013 b 

FL-M 0 S 2 

glass-supported 

-0.016 c 

-0.011 c 

Bulk MoS 2 


-0.013 a 

-0.015 a 

IL-WS 2 

glass-supported 

-0.006 d 

-0.006 d 

FL-WS 2 

glass-supported 

-0.008 c 

-0.004 c 

lL-WSe 2 

silicon-supported 

-0.0048 f 

-0.0032 f 

lL-MoSe 2 

silicon-supported 

- 

-0.0054 f 

3L-MoSe 2 

silicon-supported 

- 

-0.0045 f 


a refUH b reffl25l c refll26l d refEZI f reffl28l 


Livneh et al. reported the temperature and pressure 
dependence of the Stokes and anti-Stokes Raman spec¬ 
tra of bulk 2 H-M 0 S 2 as the energies of the A and B ex- 
citons were tuned to resonate with the exciting laserP® 
Generally, Stokes scattering peaks are stronger than anti- 


Stokes scattering peaks with a temperature-dependent 
intensity ratio, which can be exp loited practically for 
the temperature measurement! 12 ® Livneh et al. found 
that the absolute values of the temperature coefficients 
{duj/dT ) p= 0 of the frequency of nearly all of the modes 
were approximately within (1.0-1.6) x 10 -2 cm _1 /K. The 
temperature coefficient of E\ g was found to be larger 
in magnitude than that of A\ g in bulk M 0 S 2 . How¬ 
ever, in IL-M 0 S 2 , the temperature coefficient of E was 
lower in magnitude than that of A 1 , which can be at¬ 
tributed to no interlayer interactions restricting the vi¬ 
brations away from the basal plane. 12 ^Najmaei et al. re¬ 
ported laser-induced thermal effects on the Raman spec¬ 
tra of M 0 S 2 with thickness ranging from 1L to bulk. 130 
Because of the anomalous trends in frequency for E\ g 
and Ai g , they found considerable thickness-dependent 
redshifts as well as line-width changes for E\ g and A\ g 
with increasing laser power. The temperature coefficients 
of IL-WS 2 , FL-M 0 S 2 , and WS 2 in a wide temperature 
range (77-623 K) were reported by Thripuranthaka et al. 
Fig. fl5ta) and (b) show the temperature-dependent Ra- 
man spectra of the E and A 1 modes, whose frequency 
trends with temperature are summarized in Fig. pffi^c) 
and (d), respectively. 126 12 ^ Both the modes decrease in 
frequency with increasing temperature with almost the 
same temperature coefficient of about -0.006 cm -1 /K. 
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FIG. 15: Raman spectra of IL-WS 2 with the (a) E and (b) A x modes measured in a temperature range from 77 to 623 K. 
The Pos(E ) and Pos(Ai) values for IL-WS 2 are summarized in (c) and (d), respectively. Reproduced with permission from 
ref. 11271 Copyright 2014, American Chemical Society, (e) Schematic diagram of IL-M 0 S 2 on a Si 3 N 4 /Si 02 /Si substrate. Here, 
M 0 S 2 is suspended over the holes in the 20 -nm-thick SLN 4 . (f) Four typical Raman spectra of suspended IL-M 0 S 2 collected 
at 100, 180, 260, and 320 K, which are offset vertically for clarity, (g) Raman peak frequencies of the A 1 (blue squares) and 
E (red circles) modes as a function of temperature. Reproduced with permission from ref. 11251 Copyright 2014, American 
Chemical Society. 


The temperature-dependent Raman spectra of 1L- and 
FL-WSe 2 and -MoSe 2 were reported by Late et al 128 
The obtained temperature coefficients for 1L, FL, and 
bulk TMDs are summarized in Table m 

The negative temperature coefficients ( xt ) in Table 
were obtained by fitting the peak positions and tem¬ 
perature to the equations(T) = ojq + XtT, where cjo 
is the peak position at 0 K and xt is the first-order 
temperature coefficient. The intrinsic softening with in¬ 
creasing temperature can be attributed to anharmonic- 
ity of the interatomic potentials, including contribu¬ 
tions from ( 1 ) thermal expansion and ( 2 ) phonon-phonon 
coupling . 131 132 The observed linear evolution of phonon 
frequencies in IL-M 0 S 2 has been separately attributed 
to thermal expansion 125 and phonon-phonon coupling^ 24 ] 
which needs further experimental and theoretical sup¬ 
ports in a wider temperature range. Nonlinear evolu¬ 
tion of phonon frequencies in the range 70-150 K was 
detected in IL-M 0 S 2 , which was interpreted based on 
the phenomenon of optical phonon decay of two acous¬ 
tic phonons with equal energies because of lattice po¬ 
tential anharmonicityP 33 In addition, the temperature- 


dependent frequency shifts strongly depend on the cou¬ 
pling of TMDs with various substrates 

Yan et al. observed linear softening of both E' and A[ 
in suspended IL-M 0 S 2 with increasing temperature and 
laser power under 0.25 mW before heating saturation . 125 
They obtained the thermal conductivity k = 34.5 ±4 
W/mK at room temperature based on the linear temper¬ 
ature and power-dependent coefficients of the A ± mode. 
This value is less than that reported for IL-M 0 S 2 pre¬ 
pared by CVD (^52 W/mK ) 135 and much less than the 
ultrahigh thermal conductivity of graphene (4800-5600 
W/mK) Cai et al predicted a ft of 23.2 W/mK for 1 L- 
M 0 S 2 by nonequilibrium Green’s function calculations, 
which is in goo d ag reement with the reported value of 
34.5 ±4 W/mK . 137 They also calculated the Grffneissen 
parameter ( 7 ) for all of the acoustic and optical modes at 
the T point. The positive 7 value suggests a positive co¬ 
efficient of thermal expansion even at low temperature, 
which is in contrast to the negative thermal expansion 
observed in graphene at low temperatures because of the 
negative 7 of the ZA mode.EMlM] 

Furthermore, local heating in TMD flakes induced by 
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laser power has been reported in both supported and sus¬ 
pended IL-M 0 S 2 , where both E and A 1 soften with in¬ 
creasing laser power (see Fig. S5, ESI) e anc [ ^ 

linearly softened with increasing power under 0.25 mW, 
but saturated in the range 0.25-0.8 mW (see Fig. S5, 
ESI)™ The appearance of nonlinear effects is because 
of either the nonlinearity of absorption or high orders of 
the temperature-dependent coefficients. Moreover, dif¬ 
ferences in the thermal expansion coefficients between the 
TMD flake and substrate could lead to strain in TMDs, 
which further enhances the softening effect. 


C. Phonon renormalization under an electric field 

A IL-M 0 S 2 transistor can exhibit an on/off ratio of 
^10 8 and an electron mobility of ^200 cm -2 V -1 s -1 P^ 
The electronic mobility of semiconductors is limited by 
the strong electron-phonon interactions, which will also 
have a significantly affect on the phonon frequencies. Re¬ 
cently, Chakraborty et al. observed the phonon renormal¬ 
ization of in a IL-M 0 S 2 transistor. 111 The IL-M 0 S 2 
field-effect transistor (FET) and its device performance 
are shown in Fig. [l6|^a). Fig. E§b) shows the evolution 
of E and A 1 of IL-M 0 S 2 at different top-gate voltages, 
whose frequencies and FWHMs are summarized in Fig. 
16 A). At a maximum voltage of 2.0 eV, Pos(A 1 ) soft¬ 
ens by 4 cm -1 , compared with only ^0.6 cm -1 for E . 
The FWHM of A 1 significantly increases by ^6 cm -1 , 
whereas E does not significantly change. The gate in¬ 
duced electron concentration (n) can be estimated by 
ne = Ctg(Vtg~Vt ), where Vt is the threshold gate volt¬ 
age and Ctg is the top-gate capacitance. Thus, the A 1 
mode frequency shift can be used as a readout of the car¬ 
rier concentration in M 0 S 2 devicesP^Tl Chakraborty et al. 
calculated the electron-phonon coupling (EPC) based on 
first-principles density functional theory, which showed 
that A 1 couples much more strongly with electrons than 
the E' mode; 69 111 The structure distortions due to A[ 
do not break the symmetry of M 0 S 2 , which means that 
all of the electronic states can have a nonzero value in 
the electron-phonon matrix. Electron doping leads to 
occupation of the bottom of the CB at the K point, 
which results in a significant change in the EPC (A n ) 
of Ai g . In contrast, the matrix element vanishes for E , 
resulting in a weak dependence on doping. The phonon 
linewidth at a specific n (T n ) is composed of two terms: 
r n = r EPC + r an , where r EPC is the EPC contribu¬ 
tion and T an is the contribution from anharmonic ef¬ 
fects. The ratio of r n ^o/r n =o is found to follow the 
calculated strengthening EPC value A n ^o/A n =o> which 
indicates that the increase in the FWHM of the A mode 
(Fig. [Rife)) is because of the increase in the EPC value 
with doping. 111 ! 

It should be noted that A 1 renormalization is differ¬ 
ent from G-mode renormalization in graphene. The lat¬ 
ter is because of phonon-induced electron-hole (e-h) pair 


creation. 140 The phonon-induced e-h pairs are blocked 
in doped graphene, which affects the phonon self-energy. 
The G-mode renormalization in graphene is the break¬ 
down of the adiabatic approximation, 141 but A 1 renor¬ 
malization occurs within the adiabatic approximation. 


D. Substrate effect 


TMDs are usually prepared by micromechanical exfo¬ 
liation and transferred onto a SiC^/Si substrate. How¬ 
ever, much higher conducting (e.g., Au and few-layer 
graphene (FLG)) and insulating (e.g., Gel-Film®, h-BN 
flakes, and muscovite mica flakes) samples can be used 
as substrates!^ Fig. 17'a) shows E' and A[ of IL-M 0 S 2 


on six different substrates. The intensities of the E and 
A 1 modes are clearly modulated by the substrate, which 
account for optical interference and absorption effects, as 
shown in Fig. EZlM for M 0 S 2 on the Si 02 /Si substrate. 
The Raman intensity of IL-M 0 S 2 on Si02 (285 nm)/Si 
is expected to be stronger than that on both Gel-Film 
and 15 nm FLG, as shown in Fig. pT^c). It should be 
noted that this interference effect of the substrate does 
not shift the peak position. However, A 1 shows a size¬ 
able stiffening of up to ^2 cm -1 with varying substrate 
in the sequence Si02, Au, Gel-Film, FLG, mica, and 
h-BN, while E is barely affected by the substrate, as 
shown in Fig. pT^a). There are two possible explanations 
for this stiffe ning : charge transfer between IL-M 0 S 2 and 
the substrate,® 1 and the dipolar interaction between 1L- 
M 0 S 2 and the fixed charges in the different substrates PH 
The temperature behavior of A -, and E of IL-M 0 S 2 
also depends on the substrateThe temperature shift 
of A 1 is nonlinear in IL-M 0 S 2 on sapphire and Si02 (Fig. 


17 M)), while that of E is almost linear. This nonlinear¬ 
ity can be attributed to the chemical bonding between 
the film and the substrate, and the mismatch in thermal 
expansion between the film and the substrate. 


E. Pressure-induced semiconductor to metallic 
transition 

TMDs that undergo lattice distortion by applying ex¬ 
ternal pressure are theoretically predicted to undergo 
a semiconductor to semimetal transition, or even full 
metallizationfS^HSl The metallization arises from the 
overlap of the valance and conduction bands as the inter¬ 
layer spacing decreases. Recently, both X-ray diffraction 
and Raman spectroscopy of ML (at ~19 GPa) and bulk 
M 0 S 2 (between 20 and 30 GPa) under high pressure con¬ 
firmed the above transitions. 99 12 GSEI 

Fig. p~8p) shows the Raman spectra of bulk M 0 S 2 at 
various pressures up to 57 GPa in both compression (de¬ 
noted by c) and decompression (denoted by d) runs. 120 
Both E\ g and E\ g develop new split-off features with 
applied pressure above 23.0 GPa. The splitting of E\ g 
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FIG. 16: (a) Schematic diagram of the IL-M 0 S 2 FET and its device performance, (b) Raman spectra of IL-M 0 S 2 at different 
top-gate voltages Vtg- Open circles are the experimental data points, the gray (red) lines are Lorentzian fits to the total 
spectrum, and the black lines are the Lorentzian fit to individual peaks, (c) Change of phonon frequency A w (top) and FWHM 
(bottom) of A x and E as a function of Vtg- Reproduced with permission from ref. 11111 Copyright 2012, American Physical 
Society. 




Raman shift (cm -1 ) Wavelength (nm) Temperature (°C) 

FIG. 17: (a) Normalized Raman spectra of IL-M 0 S 2 on mica, FLG, Gel-Film®, Au, Si02, and h-BN, which are shifted 
vertically for clarity. The dots are the experimental points and the solid lines are Lorentzian fits. The dashed solid lines 
correspond to Pos (E ) and Pos(A 1 ) on Si02- Reproduced with permission from ref. 11421 Copyright 2014, Royal Society of 
Chemistry, (b) Model geometry of optical interference for a heterostructure based on IL-M 0 S 2 . The spacer layer can be Au, 
FLG, h-BN, or mica, (c) Enhancement factor (T -1 ) as a function of wavelength for IL-M 0 S 2 on three different substrates. 
The dashed black line indicates T _1 = 1. Temperature dependence of (d) E and (e) A ± of the bulk, CVD-grown IL-M 0 S 2 
on sapphire (lL-CVD-sa), the transferred layer of CVD-grown IL-M 0 S 2 on Si02/Si (1L-CVD-Si02), and IL-M 0 S 2 on Si02/Si 
by direct mechanical exfoliation from bulk M 0 S 2 (1L-ME-Si02). Reproduced with permission from ref. 11431 Copyright 2014, 
Royal Society of Chemistry. 


has been reported to begin at 19.1 GPaP^ Furthermore, 
E\ g is associated with interlayer stacking. Its splitting 
indicates the existence of a new phase because of layer 
sliding. This new phase at high pressure is denoted as 


2 H a: whose lattice parameter c is less than that of the 
normal 2 H phase (denoted by 2 H c in ref. !l20 ). 2 H a is a 
metallic state, as revealed by the temperature-dependent 
resistivity. The two separate peaks for both E\ g and E\ g 
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FIG. 18: (a) Raman spectra of M 0 S 2 at various pressures up to 57 GPa in both compression (denoted by c) and decompression 

(denoted by d) runs. The insets represent the in-plane modes E\ g , Ei g , and E\ g , and the out-of-plane mode A\ g . (b) Pos(Uf 5 ), 
(c) Pos(Eig), (d) Pos^E^g) and Pos(Aig), and (e) Pos(Ai^)-Pos(F^ 2 5 ) of M 0 S 2 as a function of pressure. The solid and open 
symbols are the data for the compression and decompression runs, respectively. The experimental (refs. 1991144!) and theoretical 
(ref. 145) data are shown for comparison. 2H C and 2H a corresponds to the low-pressure and high-pressure phases, respectively. 
Reproduced with permission from ref. 11201 Copyright 2014, American Physical Society. 


in bulk 2 H c completely disappear at around 40 GPa, in¬ 
dicating that the transformation from 2 H c to 2 H a is com¬ 
plete. Upon decompression, the Raman modes from 2 H 
appear, which indicates the reversibility of the 2 H c and 
2 H a phases. Fig. [I8|(b), (c), and (d) summarize the data 
of the pressure dependence of the frequencies of E\ g , Ei g , 
A \ g , and E\ g . Discontinuities occur for Pos(Ef^) and 
Pos(£ 1 2 0 ) at the beginning of the phase transition (about 
20 GPa), but not for Pos^) and Pos(i4^)™BH3H 
The discontinuities of Pos(Ef p ) and Pos^^), as well 
the frequency difference Pos(Ai^) — Pos^E^g) (Fig. 
e)), clearly divide the phase into three parts, i.e., 2 H C: 
2iJ a , and the coexistence of 2 H c and 2 H a . 


VII. APPLICATIONS OF ULF RAMAN 
SPECTROSCOPY IN 2D LAYERED MATERIALS 

A. Substrate-free characterization of layer number 
of multilayer TMDs 

The properties of 2D TMDs and other layered materi¬ 
als significantly depend on their layer number (V), which 
is important for their potential applications. A reliable 
and rapid characterization technique is required to deter¬ 
mine the layer number of ultrathin 2D TMDs. Methods 
such as optical contrast, atomic force microscopy (AFM), 
PL, and Raman spectroscopy have been widely used as 
characterization tools. 

Optical contrast is the mos t po werful characterization 
tool for 1L and ML graphene 1 — and TMDJ^. Optical 
contrast (5(A)) is defined as l-I/z a fce(A)/I su &(A), where 
lake (A) and I SU &(A) are the reflected light intensity at 


the flake and substrate, respecti vely! 48 1 15( 3 S 51 t q precisely 
identify V, the experimental optical contrast must be 
compared with the theoretical contrast for different layer 
numbers. In the calculation of optical contrast, the num¬ 
ber aperture of the used objective, the thickness of the 1L 
flake, the thickness of the dielectric layer (usually Si 02 ) 
on the substrate, and the complex refractive index of the 
flakes need to be known to obtain a reliable theoreti¬ 
cal result. However, the refractive index and monolayer 
thickness of many novel layered materials are unknown. 
The thickness of the Si02 layer provided by a company 
usually has an error of ±5%. In this case, the optical 
contrast is not very useful for layer-number identifica¬ 
tion of novel layered materials. In Fig. [l9ja), we show 
the evolution of the optical contrast of 1L to 8L M 0 WS 2 
deposited on 90-nm-thick Si02/Si. It is obviously diffi¬ 
cult to match the theoretical curve with the experimental 
curve to identify its layer number because the refractive 
index of the flakes is sensitive to their layer number. 

Atomic force microscopy is a direct and powerful tech¬ 
nique to identify the layer number with a precision of 
5%P2 However, the different interactions of the AFM 
tip with the flake and the substrate will lead to discrep¬ 
ancies in the layer number. 153 In addition, the precision 
of the AFM measurement is dependent on the substrate 
roughness and the cleanliness of the flake surface. In 
Fig. pR^b), 4L and 5L M 0 WS 2 can be clearly identified 
in the optical image. However, accidental glue spots, the 
rough surfaces, and the substrate make it difficult to ac¬ 
curately determine the layer number of this real flake in 
Fig. [T^b), as shown by its AFM image (Fig. [l9|(c)) and 
the height profile across 4L and 5L M 0 WS 2 (Fig. [l9|(d)). 

Photoluminescence is based on the evolution of the 
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FIG. 19: (a) Optical contrast as a function of N for NL-M 0 WS 2 supported on 90-nm-thick Si02/Si substrate, (b) Optical and 

(c) AFM images of 4L and 5L M 0 S 2 . (d) Height profile along the dashed green lines in (c). Raman spectra of (e) 1L-7L and 
(f) 6L-11L M 0 WS 2 . The dashed lines in (e) and (f) are guides for the eye. 


peak position or peak intensity of PL peaks with flake 
thickness!^ 9 M 0 S 2 changes from an indirect bandgap 
structure in the bulk to a direct bandgap structure in 1L- 
M 0 S 2 , which also occurs for WS 2 , WSe 2 , and MoSe 2 p^ 
This transition leads to enhanced PL emission in 1L- 
M 0 S 2 and decreasing PL intensity with increasing N. The 
PL signal is normally used for rapid characterization of 
1L- and 2L-TMDs with high crystal quality. Because PL 
emission is very sensitive to temperature, defects, and 
doping, it is not possible to accurately identify the layer 
number for thick flakes or for 2D LMs without PL emis¬ 
sion. 

Raman spectroscopy has been used as a nondestruc¬ 
tive and high-through put charact erization method for 1 L 
and ML graphene) 17 * 84 * 1Q7 ^ ^^ an d shows potential for 
general 2D layered materials:- 8 75 323 The anomalous fre¬ 
quency trends for E\ g and A\ g modes in 1 L- 6 L and bulk 
MoS 2 E3as shown in Fig. [IJ can be used for layer thick¬ 
ness determination for 1 L- 5 L M0S2 and WS2. However, 
it is not sensitive enough for M0S2/WS2 flakes with layer 
number greater than 5 . In addition, for other 1 L- and 
ML-MX2, the frequency difference between E\ g and A\ g 
needs to depend on the layer number for reconstruction, 
during which the layer number must be initially deter¬ 
mined by other methods. 

It is essential to develop a substrate-free approach to 
identify the layer number of 2D LM flakes that is not 
associated with the refractive index, monolayer thickness, 
PL emission, and high-frequency optical modes of the 
material. As discussed in Section 4.2, the peak positions 
and number of C and LB modes depend on the layer 


number of MX 2 flakesP 59 In the following, we will discuss 
how to identify the layer number of M 0 WS 2 alloys whose 
physical properties are unknown. 

Using mechanical exfoliation, we obtained a series of 
M 0 WS 2 Hakes with different thicknesses.® 0 ^ When we 
measured the ULF Raman modes of bulk M 0 WS 2 , we 
determined the C mode to be at 28.8 cm -1 , as shown in 
Fig. ll9Fe). There was no Raman signal in the C mode 
range for the flake with the smallest optical contrast. It 
was assigned as IL-M 0 WS 2 because the C mode does 
not exist in IL-M 0 WS 2 . Another flake only showed two 
ULF Raman modes at 20.4 and 35.7 cm -1 . Note that, 
the C mode at 20.4 cm -1 is exactly l/y/2 times that of 
the bulk sample (28.8 cm -1 ). According to the MCM 
in Section 4.2, the flake is 2 L-M 0 WS 2 and the mode at 
35.7 cm -1 is the LB mode of 2 L-M 0 WS 2 . Then, the 
C and LB frequencies of NL-M 0 WS 2 follow the equa¬ 
tions u(C£)(N) = C(2)y/l + cosfjr/N) (N > 2) and 
w(LB 2 -)(J\0 = LB( 2) A — cos (tt/N) (.N > 2), respec¬ 
tively. We measured the ULF Raman modes of thicker 
M 0 WS 2 flakes (N = 6 —11), as summarized in Fig. [l9je). 
By comparing the theoretical frequencies of the C and LB 
modes with the experimental frequencies, we identified 
the layer number of the thick M 0 WS 2 flakes. However, 
as shown in Fig. [l^e), when N > 7, c<j(I/ j B 2 _ )(A/') ex¬ 
hibits a different lineshape to the Lorentzian one, and 
the frequency difference of cj(C f 2 ")(A/') between sequen¬ 
tial layers will be too small to distinguish each layer be¬ 
cause of the limit of spectral resolution. The phonon 
branches of C 2 and LB^ cannot be used for layer iden¬ 
tification for N > 7. As shown in Fig. Ptf), we can 








































23 



-60 -40 -20 0 20 40 60 

Raman Shift (cm -1 ) 


FIG. 20: ULF Raman spectra of t(3+3)L, and 2H-stacked 
3L- and 6L-M0S2. The inset shows an optical image 
of t(3+3)L- and 3L-M0S2. Dash-dotted lines indicate 
Pos(LB 63 ), Pos(C 63), and Pos(LB 65 ) in 2H-stacked 6L-M0S2. 


denoted as t(3+3)L M 0 S 2 . In NL-M 0 S 2 , there are N- 
1 C modes, which we denote as C NN-i, with i= 1 , 

N-l. In this notation, for a given N, the label Cwi is 
used for the highest frequency C peak. Similar notation 
can be used for the N-l LB modes in NL-M 0 S 2 . An op¬ 
tical image of t(3+3)L is shown in the inset of Fig. 
whose 3L constituent can be verified by its C (C31 at 29.2 
cm -1 ) and LB (LB 32 at 28.1 cm -1 ) modes, as shown in 
Fig. [20j Fig. [20] also shows the Raman spectrum of 
2 H-stacked 6 L-M 0 S 2 , where C 6 i, C 63 , LB 63 , and LB 65 
are detected, in agreement with the results in Fig. [3] 
and [ 6 ] However, t( 3 + 3 )L-MoS 2 exhibits a quite different 
Raman spectrum to both 3L- and 6 L-M 0 S 2 . LB 63 and 
LB 65 are observed in t( 3 + 3 )L-MoS 2 with lower frequency 
than those in 6 L-M 0 S 2 . However, the corresponding C 6 i 
and C 63 modes in 6 L-M 0 S 2 are absent in t( 3 + 3 )L-MoS 2 
and C 62 is present. Based on the linear chain model, we 
found that the interlayer coupling constant at the twisted 
interface (cq 1 ) follows a^/<ao“=75%, where is the in¬ 
terlayer coupling constant in 2 H-stacked M 0 S 2 for the LB 
modes. This indicates that the interlayer coupling is re¬ 
duced at the twisted interface. Obviously, the LB modes 
here are themselves the direct probe of the interface cou¬ 
pling in folded (or twisted) TMDs. 
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consider other phonon branches of C§ and LBq , and 
the corresponding C and L B frequencies fo llow the equa¬ 
tions uj(Cq)(N ) = C( 2)^/l + cos(3tt/N)(N > 6 ) and 
w(LB 6 -)(JV) = LB(2)^/I — cos( 37 r/ 7 V) (N > 6 ), respec¬ 
tively. The phonon branches of C \ 0 and LB{ 0 can also 
be considered for higher layer-number identification. In 
fact, the above approach can be used for any 2D LMs 
once the C and LB modes of the material are detected. 


B. Probing interface coupling in folded TMDs 

Interface coupling is emerging as a powerful method 
to engineer the physical properties of atomically thin 2D 
materials. New phy sical phenomena, including van H ove 
singularities }^!^ 63 Fermi velocity renormalization}®^! 
and unconventional quantum Hall effectshave been 
revealed in graphene bilayers and graphene/boron ni¬ 
tride heterostructures. Several attempts has been made 
to fold M 0 S 2 bilayers to investigate how interface cou¬ 
pling will af fect its optical emission, excitons, and valley 
polarization . 166 The indirect bandgap significantly 
varies with the stacking configuration, and shows the 
largest redshift for A A- and AB-( 2 H) stacked bilay¬ 
ers, and a significantly smaller but constant redshift 
for all other twist angles !® 7 Furthermore, strong valley 
and spin polarization has been achieved in folded M 0 S 2 
bilayers P® 

From the results in Section 4 . 2 , C and LB modes are 
directly related to the interlayer coupling, which can be 
used to probe the interface coupling in folded TMDs. 
Here, we present the C and LB modes in 3L-M0S2 folded 
by another 3L-M0S2 on top. This twisted 6 L M0S2 is 


VIII. PERSPECTIVE OF THE STUDY OF 
LAYERED MATERIALS 

A. Other 2D materials beyond TMDs 

Similar to the research on graphene-like materials, 
TMDs have received increasing interest, which will 
prompt studies on other 2 D materials. Geim and Grig¬ 
orieva built a 2D library including three families: the 
graphene family (graphene, hBN, graphene oxide, etc.), 
2 D chalcogenides (M 0 S 2 , WS 2 , MoTe 2 , GaSe, Bi 2 Se 3 , 
etc.) and 2 D oxides (WO3, layered Cu oxides, Mn 02 , 
etc.), as shown in Fig. [21^ a) !®^ Although several 
achievements have been made with B^Ses/B^TesP^HES 
GaTeP® GagP®, GaS^PU, InS<™ZZI, CuStP® and 
ReS 2 ^H a large number of these 2D materials, especially 
at the ultrathin scale, are still to be investigated. 

Raman spectroscopy will give information about the 
unique properties of these 2D materials. For example, in 
contrast to TMDs, ReS 2 has been found to remain a di¬ 
rect bandgap material with decreasing thickness from the 
bulk to the monolayer, and its Raman spectrum shows 
no dependence on the layer number (Fig. |22^a)), which 
is because of electronic and vibrational decoupling . 170 
The topological insulators Bi 2 Se 3 and Bi 2 Te 3 show low- 
frequency modes that are different from the interlayer C 
and LB modes, as shown in Fig. [22^b). The out-of-plane 
A\ g mode located at ^72 cm -1 in bulk Bi 2 Se 3 exhibits 
a sensitive redshift with decreasing thickness because of 
phonon softening . 173 A significant broadening of the in¬ 
plane Eg mode is also observed with decreasing thick¬ 
ness, resulting from a decrease of the phonon lifetime 
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Layered semiconductors: 
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Perovskite-type: 

LaNb 2 0 7 , (Ca,Sr) 2 Nb 3 O 10 , 
Bi 4 Ti 3 0 12 , CagTajTiO^ and so on 
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Ti0 2 , MnO z , V 2 O s , 
Ta0 3 , Ru0 2 and so on 
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FIG. 21: (a) Current 2D library. Monolayers that are stable under ambient conditions (room temperature in air) are shaded 

blue, those that are probably stable in air are shaded green, and those that are unstable in air but may be stable in an inert 
atmosphere are shaded pink. Grey shading indicates 3-D compounds that have been successfully exfoliated down to monolayer 
thickness, (b) Building van der Waals heterostructures. Reproduced with permission from ref. 11691 Copyright 2013, Nature 
Publishing Group. 
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FIG. 22: (a) Raman spectra of bulk (blue) and 1L (red) ReS 2 . Reproduced with permission from ref. 11701 Copyright 2014, 

Nature Publishing Group, (b) Raman spectra of bulk Bi 2 Se 3 and Bi 2 Te 3 . (c) Polarized Raman spectra of bulk GaSe (top) and 
a 29-nm-thick flake (bottom). Reproduced with permission from ref. 11711 Copyright 2014, American Chemical Society. 


of the intralayer modes, possibly because of enhanced 
electron-phonon couplingFig. 22 [c) shows the Ra¬ 


man spectra of thin (29 nm) and bulk GaSe in both XX 
and YX configurations. 17 ^ In bulk GaSe (upper panel), 
E\ g and E^ g are present in both XX and YX configu¬ 
ration, while A\ g is only present in the XX configura¬ 
tion. The polarization and peak position of each mode 
in FL-GaSe (29 nm) is similar to that in bulk GaSe, in¬ 
dicating a good crystal structure. However, a significant 
decrease in intensity of Raman vibrations can be seen 
from the signal-to-noise ratio of the Raman spectra of 
FL-GaSe flakes (lower panel), which is different to the 


case of M 0 S 2 , where the Raman spectrum of 45 L-M 0 S 2 
(^29 nm) is to the same as that of the bulk. 


B. van der Waals Heterostructures 

Various 2D materials can be ” restacked” or assembled 
vertically in a chosen sequence (Fig. |2l{b)) to form var¬ 
ious hybrids and heterostructures, creating materials on 
demand, which will provide fin e-tun ed properties and 
new opportunities for devices P 79 ^ 1 8 11 These hybrids and 
heterostructures are referred to as van der Waals het- 
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erostructures, where strong covalent bonds provide in¬ 
plane stability of 2D crystals, whereas relatively weak 
van-der-Waals-like forces are sufficient to keep the stack 
together. 

The graphene/MoS 2 heterostructure firstly synthe¬ 
sized by Chang et al. 182 has stimulated investigation of 
its application in electronics and optoelectronics because 
of the semiconductivity of M 0 S 2 and the high conductiv¬ 
ity of graphene. Interesting studies have been performed 
to investigate heterostructures fabricated with various 
TMDs. Before the extensive production of TMD het- 
erostructures, several theoretical works have been per¬ 
formed to investigate their electronic structuresj 35 l 183 H IM] 
For example, Kosmider et al. found that the bandgap in 
the IL-M 0 S 2 /IL-WS 2 heterostructure is direct, which is 
in contrast to both 2 L-M 0 S 2 and 2 L-WS 2 . Addition¬ 
ally, the lowest energy electron and highest energy hole 
states in the optically active K point are localized on dif¬ 
ferent monolayers, i.e., electron-holes pairs are spatially 
separated Kang et al calculated the Moire pattern 
effect on the electronic structure of the IL-M 0 S 2 /IL- 
MoSe 2 heterostructure, and found that the VBM state 
is strongly localized while the CBM state is only weakly 
localized. 183 The number of experimental studies is grow - 
ing, and preliminary results have been obtained. 187 193 



FIG. 23: (a) Optical image of a graphene flake containing 
1LG, 3LG, and t(l+3)LG. (b) Optical contrast of t(l+3)LG 
and Bernal-stacked 4LG. (c) Raman spectra of 1LG, 3LG, 
and t(l+3)LG in the C and G peak range excited by 1.96 eV. 
The G + mode is clearly identified by the fit to the t(l+3)LG 
G band. Reproduced with permission from ref. 1811 Copyright 
2014, Nature Publishing Group. 

Interlayer coupling, strain, charge transfer, lattice vi¬ 
brations, and interface diffusion within heterostructures 
and hybrids will change the band structure, Fermi level, 
and band offsets. Raman scattering is ideal to probe 
the difference of the properties from those of the con¬ 
stituent materials. Indeed, Raman scattering has been 
applied to twisted BLG to reveal changes in the band 
structure and the appearance of angle-dependent van 
Hove singularitiesP 63 Because the interactions between 
the different layers of heterostructures and hybrids are 
weak, because van der Waals forces hold them together, 


the vibrations of heterostructures and hybrids will con¬ 
sist of those of the individual building blocks. Therefore, 
Raman spectroscopy is useful to probe the stoichiometry 
of heterostructures and hybrids. For example, the Ra¬ 
man peaks of graphene and FeCls can be identified in 
FeCl 3 -intercalated FLGs using Raman spectroscopy^^ 
while at the same time probing the layer coupl ing, sta¬ 
bility, charge transfer, strain, and orientationP^ The 
E and A 1 modes of IL-M 0 S 2 and IL-WS 2 are all ob¬ 
served in the Raman spectra of the M 0 S 2 /WS 2 het¬ 
erostructure, indicating formation of a IL-M 0 S 2 /IL-WS 2 
heterostructureP 88 The C and LB modes are mainly de¬ 
termined by the weak interlayer interactions, and will 
thus be an effective too l to probe the interlayer coupling 
in heterostructuresRaman spectroscopy can also be 
used to monitor/optimize the growth conditions for het¬ 
erostructures and hybrids. 196 


Several advances in ULF Raman spectroscopy 
have been made in investigating twisted multilayer 
graphene JHH 197 Fig. [23^ a) shows the optical image 
of a flake containing monolayer graphene (1LG), tri¬ 
layer graphene (3LG), and twisted (1+3) layer graphene 
(t(l+3)LG). Here, t(m + n)LG refers to artificially as¬ 
sembling m-layer graphene (raLG, m > 1) on n-layer 
graphene (nLG, n > 1), where mLG and nLG are twisted 
by an angle 6 t with respect to each other. 6 t between 
the layers can be used to tune the o ptical and electronic 
properties of the twisted system. 161 161 198 The optical 
contrast of t(l+3)LG shown in Fig. [23^b) is quite differ¬ 
ent to that of four-layer graphene (4LG). An additional 
feature appears at around 2.0 eV for t(l+3)LG, revealing 
that the band structure of twisted t(l+3)LG is modified 
after twisting compared with Bernal-stacked graphene 
layers. Fig. [23^c) shows the Raman spectra of 1LG, 
3LG, and t(l+3)LG in the C and G peak region. There 
is an additional G + peak in t(l+3)LG, and the G peak in 
t(1+3)LG has higher intensity, suggesting weak coupling 
between the 1LG and 3LG constituents in t(l+3)LG. The 
C mode of 3LG is observed at 37 cm -1 (C31).® Another 
C mode at 22 cm -1 ( 032 )^^ is not observed for 3LG. 
t(l+3)LG contains four layers. In Bernal-stacked 4LG, 
three C mo des ar e expected at 41 (C41), 31 (C42), and 17 
(C43) cm -1 C3HS1 However, only two C modes at 22 and 
37 cm -1 are observed, whose frequency is close to that 
of 3LG. The presence of C31 and C32 in t(l+3)LG in¬ 
dicates that the weaker interlayer coupling between the 
1LG and 3LG constituents in t(l+3)LG makes the C 
mode vibration mainly localized in the 3LG constituent. 
In addition, the C31 intensity of t(l+3)LG excited by 
1.96 eV is about 100 times higher than that of 3LG. All 
of the above results demonstrate that the peak position 
and intensity of the C modes in t(m + n)LG are suit¬ 
able to probe the interlayer coupling at the interface of 
t(m + n)LGwhich paves way for fundamental research 
of van der Waals interface coupling in 2D crystal based 
heterostructures. 
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IX. CONCLUSIONS 

Two-dimensional layered TMDs have extraordinary 
electronic, spintronic, valleytronic, and optoelectronic 
properties, and appealing applications in devices. We 
have reviewed the phonon structures of 1L, ML, and 
bulk MX 2 compounds and the capability of Raman spec¬ 
troscopy to investigate the properties of various TMDs. 
Nonresonant Raman scattering can determine the layer 
number and phonon structure at and around T points, 
probe interlayer coupling, and determine the alloy com¬ 
position and degree of various 2D crystal alloys. Res¬ 
onant Raman spectroscopy has potential to determine 
the phonon structure at the M point or along nonpla- 
nar directions and the intrinsic spin-orbit coupling in 
the valence band. Both nonresonant and resonant Ra¬ 


man spectroscopy can be used to investigate the effect of 
external perturbation, such as strain, temperature, pres¬ 
sure, charge transfer, and electric field, on TMDs. Ra¬ 
man spectroscopy is expected to be an effective tool for 
TMD-based van der Waals heterostructures and other 2D 
materials beyond TMDs. In particular, the measurement 
of C and LB modes will increase our knowledge of inter¬ 
face coupling in various van der Waals heterostructures. 

Acknowledgments 

We acknowledge support from the special funds for 
Major State Basic Research of China, contract No. 
2009CB929301, the National Natural Science Foundation 
of China, grants 11225421, 11434010 and 11474277. 


Electronic address: phtan@semi.ac.cn 

1 K. S. Novoselov and A. H. Castro Neto, Phys. Scr. 2012, 
014006 (2012). 

2 F. Bonaccorso, A. Lombardo, T. Hasan, Z. Sun, 
L. Colombo, and A. C. Ferrari, Materials Today 15 , 564 
( 2012 ). 

3 J. N. Coleman, M. Lotya, A. ONeill, S. D. Bergin, P. J. 
King, U. Khan, K. Young, A. Gaucher, S. De, R. J. Smith, 
et al., Science 331 , 568 (2011). 

4 Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, 
and M. S. Strano, Nat. Nanotechnol. 7 , 699 (2012). 

5 M. Chhowalla, H. S. Shin, G. Eda, L. J. Li, K. P. Loh, 
and H. Zhang, Nat. Chem. 5, 263 (2013). 

6 M. Xu, T. Liang, M. Shi, and H. Chen, Chemical Reviews 
113 , 3766 (2013). 

7 J. A. Wilsona and A. D. Yoffea, Adv. Phys. 18 , 193 
(1969). 

8 K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, 
Phys. Rev. Lett. 105 , 136805 (2010). 

9 A. Splendiani, L. Sun, Y. B. Zhang, T. S. Li, J. Kim, 
C. Y. Chim, G. Galli, and F. Wang, Nano Lett 10, 1271 
( 2010 ). 

10 B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 
and A. Kis, Nat. Nanotechnol. 6, 147 (2011). 

11 D. Xiao, G. B. Liu, W. X. Feng, X. D. Xu, and W. Yao, 
Phy. Rev. Lett. 108 , 196802 (2012). 

12 H. L. Zeng, J. F. Dai, W. Yao, D. Xiao, and X. D. Cui, 
Nat. Nanotechnol. 7 , 490 (2012). 

13 T. Cao, G. Wang, W. P. Han, H. Q. Ye, C. R. Zhu, J. R. 
Shi, Q. Niu, P. H. Tan, E. Wang, B. L. Liu, et al., Nat. 
Commun. 3, 887 (2012). 

14 K. F. Mak, K. L. He, J. Shan, and T. F. Heinz, Nat. 
Nanotechnol. 7 , 494 (2012). 

15 K. F. Mak, K. L. McGill, J. Park, and P. L. McEuen, 
Science 344, 1489 (2014). 

16 Y. J. Zhang, T. Oka, R. Suzuki, J. T. Ye, and Y. Iwasa, 
Science 344, 725 (2014). 

17 P. H. Tan, W. P. Han, W. J. Zhao, Z. H. Wu, K. Chang, 
H. Wang, Y. F. Wang, N. Bonini, N. Marzari, N. Pugno, 
et al., Nat. Mater. 11, 294 (2012). 

18 C. Lee, H. Yan, L. E. Brus, T. F. Heinz, J. Hone, and 
S. Ryu, ACS Nano 4, 2695 (2010). 


19 X. Zhang, W. P. Han, J. B. Wu, S. Milana, Y. Lu, Q. Q. 
Li, A. C. Ferrari, and P. H. Tan, Phys. Rev. B 87 , 115413 
(2013). 

20 Y. Y. Zhao, X. Luo, H. Li, J. Zhang, P. T. Araujo, C. K. 
Gan, J. Wu, H. Zhang, S. Y. Quek, M. S. Dresselhaus, 
et al., Nano Lett. 13 , 1007 (2013). 

21 S. Najmaei, Z. Liu, W. Zhou, X. L. Zou, G. Shi, S. D. Lei, 

B. I. Yakobson, J. C. Idrobo, P. M. Ajayan, and J. Lou, 
Nat. Mater. 12, 754 (2013). 

22 G. Eda, H. Yamaguchi, D. Voiry, T. Fujita, M. Chen, and 
M. Chhowalla, Nano Lett. 11 , 5111 (2011). 

23 S. N. Shirodkar and U. V. Waghmare, Phys. Rev. Lett. 
112, 157601 (2014). 

24 L. F. Mattheiss, Phys. Rev. B 8, 3719 (1973). 

25 S. J. Sandoval, D. Yang, R. F. Frindt, and J. C. Irwin, 
Phys. Rev. B 44, 3955 (1991). 

26 W. J. Zhao, Z. Ghorannevis, L. Q. Chu, M. L. Toh, 

C. Kloc, P. H. Tan, and G. Eda, ACS Nano 7, 791 (2013). 

27 W. C. Jin, P. C. Yeh, N. Zaki, D. T. Zhang, J. T. Sad- 
owski, A. Al-Mahboob, A. M. van der Zande, D. A. Ch- 
enet, J. I. Dadap, I. P. Herman, et al., Phys. Rev. Lett. 
Ill, 106801 (2013). 

28 W. J. Zhao, R. M. Ribeiro, M. L. Toh, A. Carvalho, 
C. Kloc, A. H. C. Neto, and G. Eda, Nano Lett. 13 , 5627 
(2013). 

29 N. Scheuschner, O. Ochedowski, A. M. Kaulitz, R. Gillen, 
M. Schleberger, and J. Maultzsch, Phys. Rev. B 89, 
125406 (2014). 

30 T. Eknapakul, P. D. C. King, M. Asakawa, P. Buaphet, 
R. H. He, S. K. Mo, H. Takagi, K. M. Shen, F. Baum- 
berger, T. Sasagawa, et al., Nano Lett. 14 , 1312 (2014). 

31 Y. Zhang, T. R. Chang, B. Zhou, Y. T. Cui, H. Yan, Z. K. 
Liu, F. Schmitt, J. Lee, R. Moore, Y. L. Chen, et al., Nat. 
Nanotechnol. 9, 111 (2014). 

32 K. Albe and A. Klein, Phys. Rev. B 66, 073413 (2002). 

33 S. Tongay, J. Zhou, C. Ataca, K. Lo, T. S. Matthews, 
J. B. Li, J. C. Grossman, and J. Q. Wu, Nano Lett. 12, 
5576 (2012). 

34 T. Cheiwchanchamnangij and W. R. L. Lambrecht, Phys. 
Rev. B 85 , 205302 (2012). 

35 H. P. Komsa and A. V. Krasheninnikov, Phys. Rev. B 86, 
241201 (2012). 



27 


36 A. Ramasubramaniam, Phys. Rev. B 86, 115409 (2012). 

37 D. Y. Qiu, F. H. da Jornada, and S. G. Louie, Phys. Rev. 
Lett. Ill, 216805 (2013). 

38 C. D. Spataru, S. Ismail-Beigi, L. X. Benedict, and S. G. 
Louie, Phys. Rev. Lett. 92, 077402 (2004). 

39 F. Wang, G. Dukovic, L. E. Brus, and T. F. Heinz, Science 
308, 838 (2005). 

40 A. Chernikov, T. Berkelbach, H. Hill, A. Rigosi, Y. Li, 
O. Aslan, D. Reichman, M. Hybertsen, and T. Heinz, 
arXiv: 1403.4270 (2014). 

41 B. Zhu, X. Chen, and X. Cui, arXiv:1403.5108 (2014). 

42 G. Wang, X. Marie, I. Gerber, T. Amand, D. Lagarde, 
L. Bouet, M. Vidal, A. Balocchi, and B. Urbaszek, 
arXiv: 1404.0056 (2014). 

43 Z. Ye, T. Cao, K. O’Brien, H. Zhu, X. Yin, Y. Wang, 
S. G. Louie, and X. Zhang, Nature 513, 214 (2014). 

44 X. B. Yin, Z. L. Ye, D. A. Chenet, Y. Ye, K. O’Brien, 
J. C. Hone, and X. Zhang, Science 344, 488 (2014). 

45 K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. 
Khotkevich, S. V. Morozov, and A. K. Geim, Proc. Natl 
Acad. Sci. 102 , 10451 (2005). 

46 V. Nicolosi, M. Chhowalla, M. G. Kanatzidis, M. S. 
Strano, and J. N. Coleman, Science 340, 1420 (2013). 

47 M. M. Benameur, B. Radisavljevic, J. S. Heron, S. Sahoo, 
H. Berger, and A. Kis, Nanotechnology 22 , 125706 (2011). 

48 H. Li, J. M. T. Wu, X. Huang, G. Lu, J. Yang, X. Lu, 
Q. H. Zhang, and H. Zhang, Acs Nano 7, 10344 (2013). 

49 R. J. Smith, P. J. King, M. Lotya, C. Wirtz, U. Khan, 
S. De, A. O’Neill, G. S. Duesberg, J. C. Grunlan, G. Mo- 
riarty, et al., Adv. Mater. 23, 3944 (2011). 

50 X. Li, W. Cai, J. An, S. Kim, J. Nah, D. Yang, R. Piner, 
A. Velamakanni, I. Jung, E. Tutuc, et al., Science 324, 
1312 (2009). 

51 Y. H. Lee, X. Q. Zhang, W. J. Zhang, M. T. Chang, C. T. 
Lin, K. D. Chang, Y. C. Yu, J. T. W. Wang, C. S. Chang, 
L. J. Li, et al., Adv. Mater. 24, 2320 (2012). 

52 K. K. Liu, W. J. Zhang, Y. H. Lee, Y. C. Lin, M. T. 
Chang, C. Su, C. S. Chang, H. Li, Y. M. Shi, H. Zhang, 
et al., Nano Lett. 12 , 1538 (2012). 

53 N. Wakabayashi, H. G. Smith, and R. M. Nicklow, Phys. 
Rev. B 12 , 659 (1975). 

54 A. Molina-Sanchez and L. Wirtz, Phys. Rev. B 84, 155413 

( 2011 ). 

55 G. Y. Zhang, G. X. Lan, and Y. F. Wang, Lattice Vibra¬ 
tion Spectroscopy, 2nd ed. (High Education Press, China, 
1991). 

56 J. L. Verble and T. J. Wieting, Phys. Rev. Lett. 25, 362 
(1970). 

57 C. Ataca, M. Topsakal, E. Akturk, and S. Ciraci, J. Phys. 
Chem. C 115, 16354 (2011). 

58 J. Ribeiro-Soares, R. M. Almeida, E. B. Barros, P. T. 
Araujo, M. S. Dresselhaus, L. G. Cancado, and A. Jorio, 
arXiv: 1407.1226 (2014). 

59 X. Luo, Y. Y. Zhao, J. Zhang, M. L. Toh, C. Kloc, Q. H. 
Xiong, and S. Y. Quek, Phys. Rev. B 88, 195313 (2013). 

60 J. L. Verble and T. J. Wieting, Phys. Rev. Lett. 25, 362 
(1970). 

61 T. J. Wieting and J. L. Verble, Phys. Rev. B 3, 4286 
(1971). 

62 S. Horzum, H. Sahin, S. Cahangirov, P. Cudazzo, A. Ru¬ 
bio, T. Serin, and F. M. Peeters, Phys. Rev. B 87, 125415 
(2013). 

63 H. Terrones, E. Del Corro, S. Feng, J. M. Poumirol, 
D. Rhodes, D. Smirnov, N. R. Pradhan, Z. Lin, M. A. T. 


Nguyen, A. L. Elias, et al., Sci. Rep. 4, 4215 (2014). 

64 Y. Q. Cai, J. H. Lan, G. Zhang, and Y. W. Zhang, Phys. 
Rev. B 89, 035438 (2014). 

65 A. Berkdemir, H. R. Gutierrez, A. R. Botello-Mendez, 
N. Perea-Lopez, A. L. Elias, C. I. Chia, B. Wang, V. H. 
Crespi, F. Lopez-Urias, J. C. Charlier, et al., Sci. Rep. 3, 
1755 (2013). 

66 C. H. Chang, X. F. Fan, S. H. Lin, and J. L. Kuo, Phys. 
Rev. B 88, 195420 (2013). 

67 L. Liang and V. Meunier, Nanoscale 6, 5394 (2014). 

68 P. H. Tan, Z. Y. Xu, X. D. Luo, W. K. Ge, Y. Zhang, 
A. Mascarenhas, H. P. Xin, and C. W. Tu, Applied 
Physics Letters 89, 101912 (2006). 

69 G. L. Frey, R. Tenne, M. J. Matthews, M. S. Dresselhaus, 
and G. Dresselhaus, Phys. Rev. B 60, 2883 (1999). 

70 R. Loudon, Adv. Phys. 50, 813 (2001). 

71 E. Smith and G. Dent, The Theory of Raman Spectroscopy 
(John Wiley V Sons, Ltd, 2005), pp. 71-92. 

72 W. Zhao, Z. Ghorannevis, K. K. Amara, J. R. Pang, 
M. Toh, X. Zhang, C. Kloc, P. H. Tan, and G. Eda, 
Nanoscale 5, 9677 (2013). 

73 X. Luo, Y. Y. Zhao, J. Zhang, Q. H. Xiong, and S. Y. 
Quek, Phys. Rev. B 88, 075320 (2013). 

74 P. Tonndorf, R. Schmidt, P. Bottger, X. Zhang, J. Borner, 
A. Liebig, M. Albrecht, C. Kloc, O. Gordan, D. R. T. 
Zahn, et al., Opt. Express 21 , 4908 (2013). 

75 S. L. Li, H. Miyazaki, H. Song, H. Kuramochi, S. Naka- 
harai, and K. Tsukagoshi, ACS Nano 6, 7381 (2012). 

76 G. Plechinger, J. Mann, E. Preciado, D. Barroso, 
A. Nguyen, J. Eroms, C. Schller, L. Bartels, and T. Korn, 
Semicond. Sci. Technol 29, 064008 (2014). 

77 G. Plechinger, S. Heydrich, J. Eroms, D. Weiss, C. Schller, 
and T. Korn, Appl. Phys. Lett. 101, 101906 (2012). 

78 H. L. Zeng, B. R. Zhu, K. Liu, J. H. Fan, X. D. Cui, and 
Q. M. Zhang, Phys. Rev. B 86, 241301 (2012). 

79 M. Boukhicha, M. Calandra, M. A. Measson, O. Lancry, 
and A. Shukla, Phys. Rev. B 87, 195316 (2013). 

80 P.-H. Tan, J.-B. Wu, W.-P. Han, W.-J. Zhao, X. Zhang, 
H. Wang, and Y.-F. Wang, Phys. Rev. B 89, 235404 
(2014). 

81 J.-B. Wu, X. Zhang, M. Ijas, W.-P. Han, X.-F. Qiao, X.-L. 
Li, D.-S. Jiang, A. C. Ferrari, and P.-H. Tan, Nat Com- 
mun 5, 5309 (2014). 

82 J. Zi, K. Zhang, and X. Xie, Phys. Rev. B 55, 9263 (1997). 

83 S. Piscanec, M. Cantoro, A. C. Ferrari, J. A. Zapien, 
Y. Lifshitz, S. T. Lee, S. Hofmann, and J. Robertson, 
Phys. Rev. B 68, 241312 (2003). 

84 P. H. Tan, Y. Deng, and Q. Zhao, Phys. Rev. B 58, 5435 
(1998). 

85 E. H. Martins Ferreira, M. V. O. Moutinho, F. Stavale, 
M. M. Lucchese, R. B. Capaz, C. A. Achete, and A. Jorio, 
Phys. Rev. B 82, 125429 (2010). 

86 W. Shi, X. Zhang, X. L. Li, Q. Q. Li, J. B. Wu, and P. H. 
Tan, Unpublished (2014). 

87 H. Richter, Z. P. Wang, and L. Ley, Solid State Commun. 
39, 625 (1981). 

88 Y. Chen, D. O. Dumcenco, Y. Zhu, X. Zhang, N. Mao, 
Q. Feng, M. Zhang, J. Zhang, P.-H. Tan, Y.-S. Huang, 
et al., Nanoscale 6, 2833 (2014). 

89 H.-P. Komsa and A. V. Krasheninnikov, J. Phys. Chem. 
Lett. 3, 3652 (2012). 

90 J. Kang, S. Tongay, J. Li, and J. Wu, J. Appl. Phys. 113, 
143703 (2013). 

91 Y. Chen, J. Xi, D. O. Dumcenco, Z. Liu, K. Suenaga, 



28 


D. Wang, Z. Shuai, Y.-S. Huang, and L. Xie, ACS Nano 
7, 4610 (2013). 

92 J. Mann, Q. Ma, P. M. Odenthal, M. Isarraraz, D. Le, 

E. Preciado, D. Barroso, K. Yamaguchi, G. V. Palacio, 

A. Nguyen, et al., Adv. Mater. 26, 1399 (2014). 

93 C. Ramkumar, K. P. Jain, and S. C. Abbi, Phys. Rev. B 
54, 7921 (1996). 

94 D. O. Dumcenco, K. Y. Chen, Y. P. Wang, Y. S. Huang, 
and K. K. Tiong, J. Alloy. Compd. 506, 940 (2010). 

95 M. Zhang, J. Wu, Y. Zhu, D. O. Dumcenco, J. Hong, 
N. Mao, S. Deng, Y. Chen, Y. Yang, C. Jin, et al., ACS 
Nano 8, 7130 (2014). 

96 K. Golasa, M. Grzeszczyk, P. Leszczynski, C. Faugeras, 
A. A. L. Nicolet, A. Wysmolek, M. Potemski, and 

A. Babinski, Appl. Phys. Lett. 104, 092106 (2014). 

97 J. M. Chen and C. S. Wang, Solid State Commun. 14, 
857 (1974). 

98 T. Sekine, K. Uchinokura, T. Nakashizu, E. Matsuura, 
and R. Yoshizaki, J. Phys. Soc. Jpn. 53, 811 (1984). 

99 T. Livneh and E. Sterer, Phys. Rev. B 81, 195209 (2010). 

100 A. M. Stacy and D. T. Hodul, Journal of Physics and 
Chemistry of Solids 46, 405 (1985). 

101 J. H. Fan, P. Gao, A. M. Zhang, B. R. Zhu, H. L. Zeng, 

X. D. Cui, R. He, and Q. M. Zhang, J. Appl. Phys. 115, 

053527 (2014). 

102 J. R. Lince, S. V. Didziulis, and J. A. Yarmoff, Phys. Rev. 
B 43, 4641 (1991). 

103 B. Windom, W. G. Sawyer, and D. Hahn, Tribology Let¬ 
ters 42, 301 (2011). 

104 L. F. Sun, J. X. Yan, D. Zhan, L. Liu, H. L. Hu, H. Li, 

B. K. Tay, J. L. Kuo, C. C. Huang, D. W. Hewak, et al., 

Phys. Rev. Lett. Ill, 126801 (2013). 

105 B. Chakraborty, H. Matte, A. K. Sood, and C. N. R. Rao, 
J. Raman Spectrosc. 44, 92 (2013). 

106 G. Salien, L. Bouet, X. Marie, G. Wang, C. R. Zhu, W. P. 
Han, Y. Lu, P. H. Tan, T. Amand, B. L. Liu, et al., Phys. 
Rev. B 86, 081301 (2012). 

107 A. C. Ferrari and D. M. Basko, Nat. Nanotechnol. 8, 235 
(2013). 

108 F. Cerdeira, E. Anastassakis, W. Kauschke, and M. Car¬ 
dona, Phys. Rev. Lett. 57, 3209 (1986). 

109 S. Bertolazzi, J. Brivio, and A. Kis, ACS Nano 5, 9703 

( 2011 ). 

110 T. Korn, S. Heydrich, M. Hirmer, J. Schmutzler, and 

C. Schuller, Appl. Phys. Lett. 99, 102109 (2011). 

111 B. Chakraborty, A. Bera, D. V. S. Muthu, S. Bhowmick, 
U. V. Waghmare, and A. K. Sood, Phys. Rev. B 85, 
161403 (2012). 

112 C. H. Chang, X. F. Fan, S. H. Lin, and J. L. Kuo, Phys. 
Rev. B 88, 195420 (2013). 

113 H. J. Conley, B. Wang, J. I. Ziegler, R. F. Haglund, S. T. 
Pantelides, and K. I. Bolotin, Nano Lett. 13, 3626 (2013). 

114 M. Buscema, M. Barkelid, V. Zwiller, H. S. J. van der 
Zant, G. A. Steele, and A. Castellanos-Gomez, Nano Lett. 
13, 358 (2013). 

115 D. Voiry, H. Yamaguchi, J. W. Li, R. Silva, D. C. B. Alves, 
T. Fujita, M. W. Chen, T. Asefa, V. B. Shenoy, G. Eda, 
et al., Nat. Mater. 12, 850 (2013). 

116 C. R. Zhu, G. Wang, B. L. Liu, X. Marie, X. F. Qiao, 
X. Zhang, X. X. Wu, H. Fan, P. H. Tan, T. Amand, et al., 
Phys. Rev. B 88, 121301 (2013). 

117 J. S. Ross, S. F. Wu, H. Y. Yu, N. J. Ghimire, A. M. 
Jones, G. Aivazian, J. Q. Yan, D. G. Mandrus, D. Xiao, 
W. Yao, et al., Nat. Commun. 4, 1474 (2013). 


118 J. S. Ross, P. Klement, A. M. Jones, N. J. Ghimire, J. Q. 
Yan, D. G. Mandrus, T. Taniguchi, K. Watanabe, K. Ki- 
tamura, W. Yao, et al., Nat. Nanotechnol. 9, 268 (2014). 

119 A. A. Mitioglu, P. Plochocka, J. N. Jadczak, W. Escoffier, 

G. L. J. A. Rikken, L. Kulyuk, and D. K. Maude, Phys. 
Rev. B 88, 245403 (2013). 

120 Z.-H. Chi, X.-M. Zhao, H. Zhang, A. F. Goncharov, S. S. 
Lobanov, T. Kagayama, M. Sakata, and X.-J. Chen, Phys. 
Rev. Lett. 113, 036802 (2014). 

121 H. Sahin, S. Tongay, S. Horzum, W. Fan, J. Zhou, J. Li, 

J. Wu, and F. M. Peeters, Phys. Rev. B 87, 165409 (2013). 

122 Y. L. Wang, C. X. Cong, C. Y. Qiu, and T. Yu, Small 9, 
2857 (2013). 

123 M. Ghorbani-Asl, N. Zibouche, M. Wahiduzzaman, A. F. 
Oliveira, A. Kuc, and T. Heine, Sci. Rep. 3, 2961 (2013). 

124 N. A. Lanzillo, A. Glen Birdwell, M. Amani, F. J. Crowne, 
P. B. Shah, S. Najmaei, Z. Liu, P. M. Ajayan, J. Lou, 
M. Dubey, et al., Appl. Phys. Lett. 103, 093102 (2013). 

125 R. Yan, J. R. Simpson, S. Bertolazzi, J. Brivio, M. Wat¬ 
son, X. Wu, A. Kis, T. Luo, A. R. Hight Walker, and 

H. G. Xing, ACS Nano 8, 986 (2014). 

126 M. Thripuranthaka, R. V. Kashid, C. S. Rout, and D. J. 
Late, Appl. Phys. Lett. 104, 081911 (2014). 

127 M. Thripuranthaka and D. J. Late, Acs Appl. Mater. In¬ 
ter. 6, 1158 (2014). 

128 D. J. Late, S. N. Shirodkar, U. V. Waghmare, V. P. 
Dravid, and C. N. R. Rao, ChemPhysChem 15, 1592 
(2014). 

129 P. H. Tan, Y. Deng, and Q. Zhao, Phys. Rev. B 58, 5435 
(1998). 

130 S. Najmaei, Z. Liu, P. M. Ajayan, and J. Lou, Appl. Phys. 
Lett. 100, 013106 (2012). 

131 J. Menendez and M. Cardona, Phys. Rev. B 29, 2051 
(1984). 

132 P. H. Tan, Y. M. Deng, Q. Zhao, and W. C. Cheng, Appl. 
Phys. Lett. 74, 1818 (1999). 

133 A. Taube, J. Judek, C. Jastrzebski, A. Duzynska, 

K. Switkowski, and M. Zdrojek, Acs Appl. Mater. Inter. 
6, 8959 (2014). 

134 L. Su, Y. Zhang, Y. Yu, and L. Cao, Nanoscale 6, 4920 
(2014). 

135 S. Sahoo, A. P. S. Gaur, M. Ahmadi, M. J.-F. Guinel, and 
R. S. Katiyar, J. Phys. Chem. C 117, 9042 (2013). 

136 A. A. Balandin, S. Ghosh, W. Bao, I. Calizo, D. Tewelde- 
brhan, F. Miao, and C. N. Lau, Nano Lett. 8, 902 (2008). 

137 Y. Cai, J. Lan, G. Zhang, and Y.-W. Zhang, Phys. Rev. 
B 89, 035438 (2014). 

138 B. D. Kong, S. Paul, M. B. Nardelli, and K. W. Kim, 
Phys. Rev. B 80, 033406 (2009). 

139 D. Yoon, Y.-W. Son, and H. Cheong, Nano Lett. 11, 3227 

( 2011 ). 

140 J. Yan, Y. Zhang, P. Kim, and A. Pinczuk, Phys. Rev. 
Lett. 98, 166802 (2007). 

141 S. Pisana, M. Lazzeri, C. Casiraghi, K. S. Novoselov, A. K. 
Geim, A. C. Ferrari, and F. Mauri, Nat. Mater. 6, 198 
(2007). 

142 M. Buscema, G. Steele, H. van der Zant, and 
A. Castellanos-Gomez, Nano Research 7, 1 (2014), ISSN 
1998-0124. 

143 L. Su, Y. Zhang, Y. Yu, and L. Cao, Nanoscale 6, 4920 
(2014), 24676020. 

144 S. Sugai and T. Ueda, Phys. Rev. B 26, 6554 (1982). 

145 L. Hromadova, R. Martonak, and E. Tosatti, Phys. Rev. 
B 87, 144105 (2013). 



29 


146 H. Guo, T. Yang, P. Tao, Y. Wang, and Z. Zhang, J. 
Appl. Phys. 113, 013709 (2013). 

147 C. Espejo, T. Rangel, A. H. Romero, X. Gonze, and G.-M. 
Rignanese, Phys. Rev. B 87, 245114 (2013). 

148 A. P. Nayak, S. Bhattacharyya, J. Zhu, J. Liu, X. Wu, 
T. Pandey, C. Jin, A. K. Singh, D. Akinwande, and J.-F. 
Lin, Nat Commun 5, 3731 (2014). 

149 P. Blake, E. W. Hill, A. H. Castro Neto, K. S. Novoselov, 
D. Jiang, R. Yang, T. J. Booth, and A. K. Geim, Appl. 
Phys. Lett. 91, 063124 (2007). 

150 Z. H. Ni, H. M. Wang, J. Kasim, H. M. Fan, T. Yu, 
Y. H. Wu, Y. P. Feng, and Z. X. Shen, Nano Lett. 7, 
2758 (2007). 

151 W. J. Zhao, P. H. Tan, J. Zhang, and J. Liu, Phys. Rev. 
B 82, 245423 (2010). 

152 K. Fukuda, K. Akatsuka, Y. Ebina, R. Ma, K. Takada, 

I. Nakai, and T. Sasaki, ACS Nano 2, 1689 (2008). 

153 P. Nemes-Incze, Z. Osvath, K. Kamaras, and L. P. Biro, 
Carbon 46, 1435 (2008). 

154 F. Tuinstra and J. L. Koenig, J. Chem. Phys. 53, 1126 
(1970). 

155 R. J. Nemanich and S. A. Solin, Phys. Rev. B 20, 392 
(1979). 

156 A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, 
M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. 
Novoselov, S. Roth, et ah, Phys. Rev. Lett. 97, 187401 
(2006). 

157 S. Reich, C. Thomsen, and J. Maultzsch, Carbon Nan¬ 
otubes: Basic Concepts and Physical Properties (Wiley- 
VCH, 2005). 

158 A. Jorio, M. S. Dresselhaus, R. Saito, and G. Dressel- 
haus, Raman Spectroscopy in Graphene Related Systems 
(Wiley-VCH, 2011). 

159 Q. Y. Zhang, Y. C. Cheng, L. Y. Gan, and U. Schwingen- 
schlogl, Phys. Rev. B 88, 245447 (2013). 

160 X. F. Qiao, X. L. Li, X. Zhang, W. Shi, and P. H. Tan, 
Unpublished (2014). 

161 G. H. Li, A. Luican, J. M. B. L. dos Santos, A. H. C. 
Neto, A. Reina, J. Kong, and E. Y. Andrei, Nat. Phys. 6, 
109 (2010). 

162 M. Yankowitz, J. Xue, D. Cormode, J. D. Sanchez- 
Yamagishi, K. Watanabe, T. Taniguchi, P. Jarillo- 
Herrero, P. Jacquod, and B. J. LeRoy, Nat Phys 8, 382 
( 2012 ). 

163 K. Kim, S. Coh, L. Z. Tan, W. Regan, J. M. Yuk, E. Chat- 
terjee, M. F. Crommie, M. L. Cohen, S. G. Louie, and 
A. Zettl, Phys. Rev. Lett. 108, 246103 (2012). 

164 G. Trambly de Laissardiere, D. Mayou, and L. Magaud, 
Nano Lett. 10, 804 (2010). 

165 K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Fal’ko, 
M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, and 
A. K. Geim, Nat. Phys. 2, 177 (2006). 

166 A. Castellanos-Gomez, H. van der Zant, and G. Steele, 
Nano Res. 7, 1 (2014). 

167 K. Liu, L. Zhang, T. Cao, C. Jin, D. Qiu, Q. Zhou, 
A. Zettl, P. Yang, S. G. Louie, and F. Wang, Nat Com¬ 
mun 5, 4966 (2014). 

168 T. Jiang, H. Liu, D. Huang, S. Zhang, Y. Li, X. Gong, Y.- 
R. Shen, W.-T. Liu, and S. Wu, Nat Nano 9, 825 (2014). 

169 A. K. Geim and I. V. Grigorieva, Nature 499, 419 (2013). 

170 S. Tongay, H. Sahin, C. Ko, A. Luce, W. Fan, K. Liu, 

J. Zhou, Y. S. Huang, C. H. Ho, J. Y. Yan, et ah, Nat. 
Commun. 5, 3252 (2014). 

171 P. Hu, J. Zhang, M. Yoon, X.-F. Qiao, X. Zhang, W. Feng, 


P. H. Tan, W. Zheng, J. Liu, X. Wang, et ah, Nano Res. 
7, 694 (2014). 

172 D. S. Kong, W. H. Dang, J. J. Cha, H. Li, S. Meister, 
H. L. Peng, Z. F. Liu, and Y. Cui, Nano Let. 10, 2245 
( 2010 ). 

173 J. Zhang, Z. Peng, A. Soni, Y. Zhao, Y. Xiong, B. Peng, 
J. Wang, M. S. Dresselhaus, and Q. Xiong, Nano Lett. 
11, 2407 (2011). 

174 F. Liu, H. Shimotani, H. Shang, T. Kanagasekaran, 

V. Zlyomi, N. Drummond, V. I. Falko, and K. Tanigaki, 
ACS Nano 8, 752 (2014). 

175 P. Hu, L. Wang, M. Yoon, J. Zhang, W. Feng, X. Wang, 
Z. Wen, J. C. Idrobo, Y. Miyamoto, D. B. Geohegan, 
et ah, Nano Lett. 13, 1649 (2013). 

176 S. R. Tamalampudi, Y.-Y. Lu, R. Kumar U., R. Sankar, 

C. -D. Liao, K. Moorthy B., C.-H. Cheng, F. C. Chou, and 
Y.-T. Chen, Nano. Lett. 14, 2800 (2014). 

177 S. D. Lei, L. H. Ge, S. Najmaei, A. George, R. Kappera, 
J. Lou, M. Chhowalla, H. Yamaguchi, G. Gupta, R. Vaj- 
tai, et ah, ACS Nano 8, 1263 (2014). 

178 X. J. Wu, X. Huang, J. Q. Liu, H. Li, J. Yang, B. Li, 

W. Huang, and H. Zhang, Angew. Chem. Int. Ed. 53, 
5083 (2014). 

179 F. Bonaccorso, P. H. Tan, and A. C. Ferrari, ACS Nano 
7, 1838 (2013). 

180 S. Z. Butler, S. M. Hollen, L. Y. Cao, Y. Cui, J. A. Gupta, 
H. R. Gutierrez, T. F. Heinz, S. S. Hong, J. X. Huang, 
A. F. Ismach, et ah, ACS Nano 7, 2898 (2013). 

181 X. Huang, C. L. Tan, Z. Y. Yin, and H. Zhang, Adv. 
Mater. 26, 2185 (2014). 

182 K. Chang and W. X. Chen, ACS Nano 5, 4720 (2011). 

183 J. Kang, J. Li, S.-S. Li, J.-B. Xia, and L.-W. Wang, Nano 
Lett. 13, 5485 (2013). 

184 J. Kang, S. Tongay, J. Zhou, J. B. Li, and J. Q. Wu, Appl. 
Phys. Lett. 102, 012111 (2013). 

185 K. Kosmider and J. Fernandez-Rossier, Phys. Rev. B 87, 
075451 (2013). 

186 H. Terrones, F. Lopez-Urias, and M. Terrones, Sci. Rep. 
3, 1549 (2013). 

187 H. Fang, C. Battaglia, C. Carraro, S. Nemsak, B. Ozdol, 
J. S. Kang, H. A. Bechtel, S. B. Desai, F. Kronast, A. A. 
Unal, et ah, 111, 6198 (2014). 

188 S. Tongay, W. Fan, J. Kang, J. Park, U. Koldemir, J. Suh, 

D. S. Narang, K. Liu, J. Ji, J. Li, et ah, Nano Lett. 14, 
3185 (2014). 

189 P. Rivera, J. R. Schaibley, A. M. Jones, J. S. Ross, S. Wu, 
G. Aivazian, P. Klement, N. J. Ghimire, J. Yan, D. G. 
Mandrus, et ah, arXiv: 1403.4985 (2014). 

190 K. Liu, L. Zhang, T. Cao, C. Jin, D. Qiu, Q. Zhou, 
A. Zettl, P. Yang, S. G. Louie, and F. Wang, 
arXiv: 1406.6487 (2014). 

191 M.-H. Chiu, C. Zhang, H. W. Shiu, C.-P. Chuu, C.-H. 
Chen, C.-Y. S. Chang, C.-H. Chen, M.-Y. Chou, C.-K. 
Shih, and L.-J. Li, arXiv:1406.5137 (2014). 

192 X. Hong, J. Kim, S.-F. Shi, Y. Zhang, C. Jin, Y. Sun, 
S. Tongay, J. Wu, Y. Zhang, and F. Wang, Nat. Nano. 9, 
682 (2014). 

193 R. Cheng, D. Li, H. Zhou, C. Wang, A. Yin, S. Jiang, 
Y. Liu, Y. Chen, Y. Huang, and X. Duan, arXiv: 1403.3447 
(2014). 

194 W. Zhao, P. H. Tan, J. Liu, and A. C. Ferrari, J. Am. 
Chem. Soc. 133, 5941 (2011). 

195 A. Castellanos-Gomez, H. van der Zant, and G. Steele, 
Nano Res. 7, 1 (2014), ISSN 1998-0124. 



196 

197 


30 


W. Dang, H. Peng, H. Li, P. Wang, and Z. Liu, Nano 2041-1723. 

Lett. 10, 2870 (2010). 198 P. Moon and M. Koshino, Phys. Rev. B 87, 205404 (2013). 

C. Cong and T. Yu, Nat. Commun. 5, 4709 (2014), ISSN 



